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SECTION  II 

Introduction  and  Summary 


Thermal  Technology  Laboratory,  Inc.  under  the 
USAF  Contract  F33615-75-C-2014  is  performing  a 
comprehensive  analytical  and  experimental  program  to 
develop  design  techniques,  experimental  high  power 
low  specific  weight  transformer  systems,  and  high 
power  lightweight  DC  to  DC  converter  system  test  beds. 

The  effort  consists  of  the  following:  (1)  the 

development  of  a 10  KVA  inverter  transformer  suitable 
for  development  elsewhere  of  experimental  converter 
circuits  and  systems,  (2)  the  refinement  of  the  real 
time  models  and  computer  design  optimization  programs 
for  inverter,  pulse,  and  sinusoidal  waveshape  transformers, 

(3)  the  development  of  two  experimental  200  KVA  inverter 
transformer  systems,  (4)  the  design,  fabrication,  and 
testing  of  two  200  KVA  converter  system  test  beds, 

(5)  the  design  and  fabrication  of  a high  voltage  non-inductive 
load  for  testing  the  converters.  These  tasks  are 
summarized  below. 

In  Task  1,  a ten  (10)  KVA  experimental  inverter 
transformer  has  been  designed  to  Government  furnished 
requirements,  assembled,  and  subjected  to  preliminary 
tests.  A second  set  of  windings  have  been  prc.ided  for 
experimentation  and  the  computer  program  is  being  modified, 
to  correlate  with  the  measured  electrical  performance 
of  the  transformer. 

In  Task  2,  the  inverter,  pulse,  and  real  time  models 
of  the  transformer  are  being  refined  by  generalizing  the 
steady  state  and  real  time  models-  determining  and 
implementing  a program  procedure  which  will  accomodate 
both  models;  developing  an  improved  computer  aided  system 
for  design  optimization,  and  programing  on  site  the 
WPAFB  CDC  6600  computer'. 

In  Task  3,  two  200  KVA  inverter  transformer  systems 
will  be  designed,  fabricated  and  tested  in  conformance 
with  requirements  provided  by  the  Government.  Materials 
compatability  tests  will  be  conducted  and  special  cooling 
systems  designed  to  match  the  transformer  parameters. 

One  transformer  will  be  designed  to  work  with  SCR  type 
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inverters  and  the  other  designed  to  work  with  plasma 
switch  inverters. 


In  Task  4,  two  complete  operational  200  KVA  converter 
test  beds  will  be  designed,  fabricated,  and  tested  in 
accordance  with  the  requirements  to  be  provided  by  the 
Government.  One  converter  will  utilize  SCR ' s and  the 
other  plasma  switches.  A packaging  philosophy  will  be 
developed  and  implemented  in  the  electrical,  thermal, 
and  mechanical  designs.  Emphasis  will  be  placed  on 
safety,  minimum  weight,  thermal  and  electrical  performance. 

A detailed  test  plan  will  be  submitted  prior  to  testing. 

In  Task  5,  a noninductive  load  for  testing  the 
converters  will  be  designed,  fabricated,  and  tested.  The 
load  will  be  designed  to  minimize  changes  in  performance 
due  to  varying  duty  cycles  and  input  powers.  Safety,  and 
electrical  and  thermal  performance  will  be  emphasized. 

Summary 

Several  major  developments  have  occurred  in  this 
program.  As  is  often  the  case  in  research,  they  did 
not  occur  exactly  in  accordance  with  the  original  plan. 

The  total  program  however  is  essentially  on  schedule. 

New  techniques  for  fabrication  of  pie  wound 
transformers  were  developed  which  yielded  superior  designs. 
This  development  was  followed  by  the  development  of  computer 
aided  design  programs  for  pie  wound  transformers.  Several 
10KW  transformers  have  been  fabricated  and  subjected  to  a 
variety  of  tests.  Based  on  the  results  of  these  tests  both 
the  10KW  and  200  KW  transformers  will  be  pie  rather  than 
layer  wound.  A 10KW  transformer  rectifier  system  is 
presently  being  integrated  with  a breadboard  inverter  for 
final  verification. 

During  the  first  half  of  this  program,  a numerical 
method  was  developed  for  the  solution  of  the  nonlinear 
lumped  parameter  transformer  model.  This  model  was 
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developed  under  the  previous  contract  but  no  stable 
solution  had  been  found.  The  present  solution  is  for  a 
resistive  load  and  work  is  continuing  to  include  leakage 
inductance,  shunt  capacitance,  and  an  arbitrary  load 
impedance.  The  present  program  has  been  implemented  on 
the  HP  9830  machine. 

As  mentioned  above,  pie  wound  transformer  design 
programs  have  been  written.  These  are  for  interactive 
design  on  the  HP  9830  and  are  not  self  optimizing  routines. 

A decision  should  be  made  as  to  the  ultimate  application 
of  these  programs  before  a decision  whether  or  not  to 
implement  them  on  the  CDC  6600. 

The  materials  studies  support  of  the  200KW  transformer 
task  are  completed  with  exception  of  evaluation  of  some 
of  the  newer  fluorine  based  dielectric  fluids.  Preliminary 
designs  indicate  specific  weights  in  the  neighborhood  of 
0.05  lb/KVA  for  these  transformers.  Consideration  is 
presently  being  given  to  increasing  the  operating  temperature 
in  order  to  reduce  heat  exchanger  weight. 

In  support  of  the  general  materials  studies,  several 
improvements  were  made  in  the  magnetic  core  model  which 
permit  accurate  simulation  of  hysteresis  and  I-V  curves 
for  most  core  materials.  The  model  has  been  compared  to 
existing  models  and  found  to  be  superior  in  both  simplicity 
and  accuracy  in  all  cases. 

Computer  aided  design  of  transformers  is  facilitated 
through  two  transformer  design  optimization  programs, 

TRANSOP  and  TDOP2.  Additional  work  is  planned  to  improve 
their  versatility  and  usefulness  by  including  some  of  the 
newer  transformer  configurations,  i.e.,  pie  windings  and 
magnetic  shunts.  Improvements  in  ease  of  utilization, 
computational  efficiency  and  documentation  are  also  planned. 
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SECTION  III 
Technical  Discussion 

3.1  Summary  of  capabilities  of  original  program  (TRANSOP) 

| 

The  original  computer  program  developed  under  Contract 
F3 3615 -72 -C -1 944  for  computer  aided  transformer  design 
optimization  was  documented  in  report  No.  AFAPL-TR  75-15, 

Volume  II,  'Final  Technical  Report  on  Development  of 
Lightweight  Transformers  for  Airborne  High  Power  Supplies, 

Computer  Users  Manual'.  The  program  consisted  of  a 
j procedure,  FUNCTION  TDP,  and  its  associated  subroutines 

and  functions,  which  determined  a consistent  set  of 
dependent  design  parameter  values  for  a given  set  of 

independent  design  parameters.  In  addition,  the  program  I 

contained  a driver  procedure,  PROGRAM  TDPLOT,  to  aid  in 

optimizing  transformer  designs.  This  procedure  read 

control  and  independent  design  parameter  value  data  from 

an  appropriate  input  file  and  generated  corresponding 

outputs.  The  outputs  are  in  the  form  of  3 dimensional 

line  printer  and/or  Cal-Comp  plots  of  the  user  defined 

objective  function  versus  two  of  the  independent  design 

parameters  being  optimized.  Also,  a listing  of  the  final 

(optimized)  transformer  design  parameter  values  are 

computed  by  FUNCTION  TDP.  A procedure,  PROGRAM  TERMOP , 

was  developed  to  facilitate  utilization  of  PROGRAM  TDPLOT 

from  a remote  time  share  terminal.  This  program  prompts  5 

for  appropriate  free  formatted  control  and  independent 

design  parameter  data  input  and  generated  the  necessary 

control  and  input  data  card  images  on  a file  for  submission 

to  the  REMOTE  INPUT  FILE  as  an  INTERCOM  BATCH  job.  Output 

resulting  from  execution  of  PROGRAM  TDPLOT  as  an  INTERCOM 

BATCH  job  can  be  listed  at  a remote  time  share  terminal  or  i 

ROUTED  to  a batch  site  printer. 

FUNCTION  TDP  designs  3-phase  layer  wound,  vaporization 
cooled  rectangular-crossection,  C-core  transformers.  The 
user  must  specify:  (1)  the  conductor  cooling  rates  or 
standard  wire  gauge  (AWG),  (2)  core  lamination  width, 

(3)  core  lamination  thickness,  (4)  core  lamination 
insulation  thickness,  (5)  number  of  inside  core  laminations, 

(6)  number  of  outside  core  laminations,  (7)  insulation  ] 

dielectric  strength,  (8)  delta  or  wye  connected  primary 
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and  secondary  windings,  (9)  the  number  of  inside 
secondary  winding  layers,  (10)  coolant  boiling 
temperature,  (11)  coolant  vaporization  cooling 
coefficient  and  exponent,  (12)  conductor  insulation 
thicknesses,  (13)  insulation  thermal  resistivity 
coefficients,  (14)  resistivity  temperature  coefficients 
and  zero  degree  centigrade  resistivities  (15)  core  spacing 
factor,  (17)  densities,  (18)  primary  and  secondary, 
line-to-line  full  load  emf,  (19)  operating  frequency, 

(20)  coolant  heat  of  vaporization,  (21)  hysteresis  dissipation 
coefficient  and  exponent  for  the  Steinmetz  relationship, 

(22)  maximum  number  of  convergence  iterations,  (23)  full-load 
output  power,  (24)  eddy-current  dissipation  coefficient, 

(25)  core  saturation  flux  density,  (26)  Cooling  spacing 
factors,  (27)  minimum  spacings,  (28)  enclosure  thickness, 

(29)  enclosure  cooling  rate,  and  (30)  whether  conductors 
are  cooled  on  both  or  one  side.  The  function  initially 
assumes  that  the  transformer  efficiency  is  100%  and 
iteratively  converges  on  the  actual  efficiency  to  within 
a specified  error.  The  function  then  computes  the 
remaining  300  transformer  design  parameters. 

The  major  advantage  of  the  TRANSOP  program  is  its 
ability  to  provide  visual  information  on  parametric 
relationships  in  the  form  of  plots.  Another  useful 
capability  is  its  ability  to  consider  non  standard 
transformer  designs  having  different  size  inner  and  outer  core 
leg  cross-sectional  areas  and  different  inner  and  outer 
core  leg  winding  wire  sizes.  The  large  number  of  parameters 
make  it  possible  to  consider  and  analyze  many  other 
transformer  design  parameter  characteristics . 

Unfortunately,  such  detail  makes  the  program  quite 
large  with  a tendency  to  be  time  consuming  and  laborious 
to  utilize.  This  subsequently  lead  to  the  development  of 
much  simpler  and  easier  to  use  interactive  programs  in 
HP-BASIC  to  run  on  an  HP-9830  A programmable  calculator. 
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3.1.2  New  Transformer  Design  Optimization, 

Program  II  (TD0P2) 

A second  transformer  design  optimization  program 
(TDOP2)  has  been  implemented  on  CDC-6000  series  computing 
systems  in  Extended  FORTRAN.  This  is  a translated  and 
highly  modified  version  of  a simpler  and  correspondingly 
less  versatile  program  previously  implemented  in  BASIC  on 
an  HP-9830.  TDOP2  allows  ranging  of  specified  transformer 
parameters  over  specified  ranges  (e.g.,  core  area,  conductor 
size,  number  of  primary  and  secondary  layers,  etc.)  For 
each  combination  of  ranged  parameter  values  (along  with 
the  fixed  input  parameters) , the  program  computes  the 
design  properties  of  a transformer  by  means  of  an  iterative 
convergent  efficiency  computation.  In  addition,  bounds 
may  be  set  on  derived  transformer  parameters  (e.g., 
minimum  efficiency,  maximum  per  unit  resistance  and 
reactance,  one  or  more  physical  dimensions,  etc.). 

3. 1.2.1  Transformer  Design  Capabilities 

Transformer  design  parametric  relationships  incorporated 
in  to  TD0P2  include: 

Leakage  Inductance 

A standard  approximation  for  calculating  leakage 
inductance  L of  a transformer,  referenced  to  a secondary 
winding  of  N turns  is: 

L - 32McNZjT '/nZ  /o'9  (md  +t/3)  henries  (1) 


where 

N = Number  of  turns  in  secondary  winding 
Mc  = Mean  circumference  of  primary  + secondary  winding  (inches) 
1 = Height  of  cylinder  type  coil  (inches) 

n = Number  of  interleavings  of  primary  and  secondary  layers 
d = Distance  between  primary  and  secondary  layers  (inches) 
t = Primary  + secondary  winding  thickness  (inches) 

This  formula  has  been  verified  for  accuracy  in  previous 
transformer  designs,  and  is  used  in  the  program. 
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The  corresponding  leakage  reactance  in  ohms,  Xl,  is 
given  by: 

XI  = 2TTfL  (2) 

where  f = operating  frequency.  Hertz 

The  leakage  reactance  on  a per  unit  basis,  XL  is 
given  by: 

X/  = 3XlIsZ/Po  (3) 

where  Is  = secondary  winding  current,  ohms 
PQ  = output  power,  watts 

In  general,  the  leakage  reactance  of  most  designs 
could  be  maintained  at  about  0.02-0.03  by  maximizing  the 
number  of  interleavings,  n . 

Per  Unit  Resistance 

The  program  computes  the  total  dissipation  (copper  loss) 
Pc  in  watts  per  coil  by 

R:  = T.p  Rp  ~f~ -Z"s  R S (4) 

Where  Ip  = primary  winding  current,  amperes 
Rp  = primary  resistance,  ohms 
Is  = secondary  winding  current,  amperes 
Rg  = secondary  resistance,  ohms 

The  computed  resistance  values  include  both  AC 
resistivity  and  temperature  coefficient  of  resistivity 
effects . 

The  per  unit  resistance,  R,  is  given  by 

R-3PC/P0  (5) 
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In  general,  as  the  unit  surface  cooling  rate 
increases,  the  per  unit  resistance  increases.  A 
limitation  on  maximum  per  unit  resistance  can  prevent 
utilization  of  very  high  unit  surface  cooling  rates. 

This  factor  can  lead  to  the  ultimate  selection  of 
forced  oil  cooling  over  vaporization  cooling, when 
very  high  unit  surface  cooling  rates  are  not  required. 

Primary -to -secondary  interleaving 

The  primary  effect  of  an  increased  number  of 
primary  to  secondary  interleavings  is  to  reduce  the 
leakage  inductance.  The  program  assumes  that  for  a 
given  number  of  primary  and  secondary  layers,  the 
coils  will  be  arranged  so  as  to  maximize  the  number  of 
interleavings . 

AC  Resistivity  Effects 

As  conductor  diameter  and  operating  frequency  increase, 
skin  effects  become  more  pronounced.  The  net  result  is 
that  for  a given  conductor  diameter,  the  apparent  resistance 
of  the  conductor  is  greater  at  high  frequencies  than  the 
D.C.  resistance.  The  published  tabulated  data  for  this 
AC  resistivity  increase  were  fitted  to  the  equation 

Rac/ Roc  - 0.9908  + 0./768  [(  X -2.//)  + sf^L-Z.ll)7- + 0.283  ] (6) 

where 

RAc/RDC  = resistance  ratio,  AC  to  DC 

x = 0.271767  d sff1  (for  copper  conductors  only) 
d = conductor  diameter,  inches 
f = operating  frequency,  hertz 

The  relationship  is  theoretically  valid  only  for 
isolated  conductors,  but  it  has  been  applied  to  the  side- 
by-side  conductor  situation  existing  in  a transformer  coil. 
The  effect  has  been  included  in  the  program. 
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Duty  Cycle 

Duty  cycle  has  a significant  effect  on  the  choice 
of  cooling  system  and  the  resultant  overall  system  size 
and  weight.  For  very  low  duty  cycle  operation  an 
adiabatic  system  (no  explicit’  cooling  system)  can  be 
constructed  with  adequate  thermal  capacity  and  with 
reasonably  sized  cores  and  conductors.  Because  there 
are  no  external  cooling  system  components,  this  system 
can  have  decided  size  and  weight  advantages.  Such  a 
system,  however,  may  not  meet  higher  duty  cycle  requirements 
and  active  cooling  systems  would  have  to  be  employed.  The 
effects  of  thermal  capacity  on  necessary  cooling  rates  is 
considered  in  the  program  for  any  non  continuous  duty 
transformer. 

Load  Voltage 

The  necessary  full  load  secondary  rms  line-to-line 
emf  is  equal  to  0.74  for  a Three-Phase  Bridge  or  1.11  for 
a Single-Phase  Bridge  times  the  average  bridge  rectifier 
dc  output  emf  divided  by  the  cosine  of  the  value  of  .694 
times  the  commutation  overlap  angle,  i.e. 


Esl  - 0.14  * £0 / cos  (0.694  * oc  ) 3- PHASE 

£si  =/.//*  £o /cos  (0.(694-  * u)  /-PHASE 


where 

ESl  and  Eq  are  in  equal  units  and  is  in  radians. 

The  full  load  secondary  rms  winding  emf  is  equal  to 
the  line-to-line  emf  for  single-phase  and  3-phase  delta 
connected  secondary  and  to  the  line-to-line  emf  divided 
by  the  square  root  of  3 for  a y connected  secondary,  i.e. 


(8a) 


sw 


: A 4 S/N6LE  PHASE 


£»,  = E.u  / •■  y 


(8b) 


Load  Current 

The  secondary  rms  line  current  is  equal  to  0.816 
for  a Three-Phase  Bridge  or  1.00  for  a Single-Phase 
Bridge  times  the  average  bridge  rectified  dc  output 
current  times  the  rms  load  current  waveform  factor,  i.e, 


TSl  = 0.8/6  *Io*  Fl 


3 -PHASE 


(9a) 


zSL  - /.oo  *j:0  * fl 


/-PHASE 


(9b) 


where 

*SL  ant^  I0  are  ectual  units  and 

= 1.11  for  half  wine  wave  load  current  waveforms 
= 1.04  for  semi  square  wave  load  current  waveforms 
= 1.00  for  a constant  DC  load 

The  rms  secondary  winding  current  is  equal  to  the 
rms  line  current  for  Single-Phase  and  3-Phase  wYe  connected 
secondary  and  to  the  rms  line  current  divided  by  the  square 
root  of  3 for  a delta  connected  secondary,  i.e. 


I 
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SL 
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( 10a) 


X 


-W 


= Tsl/sTz  : A 


(10b) 
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3- 1.2. 2 Surface  Cooling  Rates  a Cooling  System 
Descriptions 

Attainable  surface  cooling  rates  (cores,  conductors, 
rectifiers  plates)  are  a function  of  the  cooling  system 
under  consideration. 

Three  basic  methods  of  cooling  can  be  specified. 

Natural  convection  cooling  relationships  are  presently 
not  included  in  the  program.  The  rate  of  heat  removal 
by  liquid  natural  convection  is  very  low  (on  the  order  of 
0.1  watts/ ( sq. in-C) . A brief  description  of  the  physical 
aspects  of  each  system  follows. 

Vaporization  Cooling 

In  a vaporization  cooled  transformer,  the  coil 
assemblies  are  formed  in  concentric  layers,  with  annular 
cooling  passaqes  between  layers.  When  immersed  in  a low 
boiling  point  coolant,  such  as  Freon  113  or  3 M's  FC  78, 
the  Joulean  heat  raises  the  coil  temperature,  the  coolant 
at  the  coil  surface  boils,  and  the  vapor  escapes  through 
the  cooling  passage  by  natural  bouyaney  forces.  The  heat 
of  vaporization  of  the  fluid,  combined  with  the  agitation 
resulting  from  vapor  bubble  flow,  results  in  the  capability 
of  extremely  high  heat  transfer  rates  at  the  conductor 
surface.  The  generated  vapor  must  either  be  expended  or 
condensed  and  returned  to  the  system.  Cost,  logistic,  and 
environmental  conditions  must  be  considered.  The  volume 
and  weight  of  a condenser  whose  capacity  matched  the  on 
time  dissipation  rate  would  generally  be  large,  even  if 
thermal  capacity  is  included.  The  use  of  a smaller 
capacity  condenser  (one  whose  capacity  equals  or  exceeds 
the  dissipation  rate  averaged  over  an  entire  cycle, 
including  off  time)  requires  a vapor  storage  facility. 

An  accumulator  could  be  used  for  this  purpose,  maintaining 
an  essentially  constant  pressure  system,  but  the  vapor 
generation  rate  generally  requires  too  great  a volume 
to  be  practical.  For  example,  a 3 MW  transformer  dissipating 
45  KW  (98.5%  efficient)  using  an  F113  vaporization  cooling 
system  would  generate  approximately  22  cubic  feet  of  vapor 
during  each  15  seconds  of  operation,  ignoring  thermal 
capacity  effects  and  heat  transfer  to  the  environment. 


The  alternative  to  a constant  pressure  accumulator 
coolinq  system  is  a constant  volume  system,  wherein  the 
s/stem  vapor  pressure  is  allowed  to  vary  with  the  fluid 
temperature;  that  is,  the  system  pressure  is  always  the 
vapor  saturation  pressure  corresponding  to  the  fluid 
temperature,  and  incipient  boiling  conditions  always 
prevail.  The  physical  housing  of  such  a system  including 
the  case  and  all  case  seals,  penetrations,  and  gasketed 
joints,  must  be  designed  to  withstand  the  maximum  internal 
pressure  an ticipated , as  well  as  the  possibility  of  negative 
internal  pressure  during  low  temperature  non-operating 
conditions . 

Vaporization  cooling  is  capable  of  unit  surface  cooling 
rates  of  1.0  to  50.0  watts/(sq  in-C) . Use  of  the  higher 
end  of  the  range  of  cooling  rates  available  with  this 
system  may  lead  to  system  efficiencies  below  the  minimum 
acceptable  value  for  some  applications.  Consequently, 
full  use  may  not  be  made  of  this  cooling  capability  to 
achieve  minimum  size  and  weight. 

Adiabatic  System 

An  adiabatic  system  relies  solely  on  the  thermal 
capacitance  of  the  system  elements  to  absorb  the  dissipated 
heat  without  excessive  temperature  rise.  No  explicit  heat 
transfer  to  the  environment  is  included;  rather,  it  is 
assumed  that  the  off  time  between  operating  intervals  is 
long  enough  to  ensure  stabilization  of  the  equipment  at 
or  near  ambient  temperatures  prior  to  the  next  operating 
interval,  by  natural  modes  of  conduction,  free  convection, 
and  radiation. 

The  obvious  advantage  of  the  adiabatic  system  is  the 
weight  and  volume  saving  associated  with  the  absence  of 
any  explicit  cooling  system.  Conductor  sizes  in  adiabatic 
systems  tend  to  be  larger  than  those  in  systems  employing 
explicit  cooling  means,  the  overall  weight  of  the  transformer 
is  greater,  but  system  weight  can  be  less,  since  a cooling 
system  is  not  required.  Rectifiers  in  an  adiabatic  system 
must  be  provided  with  a heat  absorbing  sink  in  the  form  of 
a metal  plate  of  sufficient  thermal  capacity. 


Surface  cooling  rates  are  so  low  that  for  practical 
purposes,  they  may  be  assumed  to  be  zero.  Dissipative 
elements  must  be  designed  to  have  sufficient  thermal 
capacity  to  absorb  the  generated  heat  during  an  operating 
interval  with  an  acceptable  temperature  rise.  These 
systems  are  generally  suitable  for  total  on  time  within 
an  operating  interval  of  about  60  seconds  or  less.  The 
time  between  operating  intervals  must  be  on  the  order  of 
24  hours  or  more  to  allow  a sufficient  cooling  down  period. 

Forced  Liquid  Cooled  Systems 

Forced  oil  cooling  systems  provide  a compromise 
between  vaporization  cooled  systems  and  adiabatic  systems. 
Vaporization  cooled  systems  provide  a very  high  heat 
transfer  capability,  but  require  a cooling  system 
carefully  designed  for  each  application.  Adiabatic 
systems  provide  a very  low  thermal  capability,  but  do  so 
with  extreme  physical  simplicity.  Forced  oil  cooled 
systems  provide  a moderate  thermal  capability,  using 
cooling  system  components  which  are  usually  somewhat  smaller 
than  for  vaporization  cooling  but  also  require  a pump  and 
motor.  The  tradeoffs  must  be  evaluated. 

A forced  oil  system  is  physically  similar  in 
construction  to  a vaporization  cooled  system,  although 
actual  dimensions  may  differ.  Coil  assemblies  are  formed 
in  concentric  layers,  with  annular  cooling  spaces  between 
layers.  The  coolant  is  an  oil  chosen  on  the  basis  of 
thermodynamic  properties,  dielectric  breakdown  characteristic, 
cost  and  availability,  and  other  pertinent  factors.  The 
coolant  is  forced  through  the  cooling  spaces  by  an  external 
pump.  The  transformer  design  must  minimize  flow  paths 
in  parallel  with  the  cooling  spaces  and  must  ensure  an 
adequate  oil  flow  velocity  in  each  channel  and  in  all 
portions  of  a given  channel. 

Forced  liquid  cooling  is  capable  of  unit  surface 
cooling  rates  of  0.1  to  1.0  watts/(sq  in  - C) , dependent  on 
the-  thermodynamic  properties  of  the  fluid  coolant  and  on 
the  coolant  velocity.  Coolant  flow  rate  must  be  correlated 
with  the  head-flow  characteristics  of  the  available  pumps, 
and  with  total  system  pressure  drop,  including  external 


heat  exchanger.  The  available  cooling  rates  for  moderate 
flow  rates  of  a typical  transformer  oil  can  be  compatable 
with  the  required  system  characteristics  and  minimum 
weight  and  volume. 

Forced  Air  Cooling 

Forced  liquid  cooling  relationships  are  also 
characteristic  of  a forced  gas  cooled  system  when  gas 
compression  is  not  significant. 

Forced  air  cooling,  using  ambient  air  blown  through  an 
open  coil  transformer  structure  and  an  open  rectifier  bank, 
can  achieve  cooling  rates  in  excess  of  0.1  watts/(sq  in-C) , 
and  requires  only  a blower  and  filter  as  an  external 
cooling  system.  When  very  high  voltages  are  involved, 
this  type  cooling  system  may  not  be  acceptable  on  the 
basis  of  safety  and  contamination  with  subsequent  breakdown. 

3.1.3  Transformer  Design  Optimization  Function 
Considerations 

The  obvious  and  usually  specified  objective  function 
for  a design  optimization  is  the  transformer  weight.  This 
specification  usually  arises  because  the  design  optimization 
of  the  overall  system  is  broken  into  sections  which  include 
the  design  optimization  of  subsystems  with  the  transformer 
being  treated  as  an  individual  component  whose  weight  adds 
to  the  total  weight  of  the  system.  Therefore,  by 
minimizing  the  weight  of  each  subsystem  component  the 
overall  system  weight  would  be  minimized.  It  is  believed 
that  the  preceeding  is  not  necessarily  a valid  assumption. 

It  appears  that  if  a minimum  weight  optimization  of  the 
overall  system  design  is  carried  out  on  the  overall 
integrated  system,  the  resulting  optimal  transformer  design 
may  not  be  a minimum  weight  transformer.  Therefore, 
transformer  weight  is  not  necessarily  a suitable  transformer 
design  optimization  objective  function  in  itself. 


systems  designed  to  match  the  transformer  parameters. 
One  transformer  will  be  designed  to  work  with  SCR  type 
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To  support  the  preceding  statements  and  derive  a 
suitable  objective  function  for  minimizing  the  overall 
system  weight,  the  following  general  analysis  of  a 
representative  system  is  presented. 


Fig.  3-1 

Minimum  Weight  System  Block  Diagram 

Fig.  31  is  a block  diagram  representation  of  a 
sub  system  interconnected  to  the  remainder  of  a system 
whose  overall  weight  is  to  be  minimized,  where: 

Pt  = input  power  to  sub  system  « 

P0  = useful  output  power  from  sub-system 
PL  = sub  system  power  losses 

Wy  = weight 

Cy  = specific  weight 

^y  = efficiency 

Xs  = source  . J 

Xc  = sub  system  (component) 

Xl  = load 
Xjj  = heat  sink 

The  total  system  weight  WT  is  equal  to  the  sum  of  the 
weights  of  the  sub  systems,  i.e. 

WT  = Ws  + Wc  + W4  + VJH  (11) 


15 


2 


Substituting  for  the  sub-system  weights  their  corresponding 
specific  weight  and  power,  i.e. 

W - C * P <12> 

We  can  rewrite  the  system  weight  relationship  in  terms 
of  specific  weights,  i.e. 

VJr  = £s  * Px  + Cc  * Po  + CL  * Po  + CH  * Pl  (13) 

Rewriting  the  input  and  loss  power  in  terms  of  the  output 
power  and  sub  system  efficiency,  i.e. 

pz  --  p°lnc  (14) 

PL--Po(i-')  <15> 

Substituting,  the  total  system  weight  equation  now 
becomes 

WV  =Cs*Po/nc  + Cc  *Po  +CL  *Po  + CH*P0  * (-/  + £)  (16) 

Dividing  through  by  the  sub-system  output  power  we  have 
a resulting  relationship  for  the  specific  weight  of  the 
total  system,  i.e. 

cT  = ^ = yf~c  + c<  + ^ + c»  *(ic  (17) 

The  individual  specific  weights  can  be  represented  as 
ratios  of  the  specific  weight  of  the  sub-system  of 
interest,  i.e. 

Cs  = Rs  * Cc  (18a) 

Cl  r Pl  * Cc  (18b) 

Ch  - * Cc  <18^) 


The  specific  weight  of  the  total  system  can  now 
be  specified  in  terms  of  the  specific  weight  and 
efficiency  of  the  sub-system  of  interest  and  the  ratio 
of  all  other  specific  weights  to  it  i.e., 

CT  = ^ + Cc  + rl  cc  + rh  Cc  (-/  - nc ) (19) 

Factoring  out  Cc: 


CT  ' Cc  * ( + Rl  Rh  + + 1 ) (20a) 

= ^ ( R%  + Rh  + nC  (l  + Rl-Rh)) 


This  relationship  can  be  used  as  the  objective  function 
for  the  design  optimization  of  a sub-system  (component) 
in  a minimum  weight  system.  The  objective  function  can 
be  simplified  if  certain  relationships  exist. 


IF  Rs  + Rh  ~ O OR  ^ + 1 + 

THEN  CT  <x  Cc 

Cc 

or  IF  /?s  +/?„»!+  Rl  ' Rh  THEN  C7  °< 


(21a) 


(21b) 


When  little  or  no  information  is  available  on  the  specific 
weight  ratios  it  is  felt  that  the  later  simplified  objective 
function  gives  a more  suitable  optimized  sub-system  design 
than  the  former  simplified  objective  function. 


Therefore,  in  the  design  optimization  of  a transformer 
or  transformer-rectifier  system  as  a sub-system  in  a 
system  whose  total  weight  is  to  be  minimized,  it  is  more 
desirable  to  use  weight  divided  by  efficiency  rather  than 
just  weight  as  the  design  optimization  objective  function. 
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In  general  the  proceeding  statement  says  that  it  is 
desirable  to  maximize  efficiency  at  the  same  time  weight 
is  being  minimized.  In  maximizing  efficiency,  losses  are 
being  minimized,  and  therefore  we  should  be  able  to  restate 
our  objective  as  the  desire  to  minimize  losses  and  weight. 

The  two  relationships  are  not  identical  as  they  possess 
different  weightings  as  is  indicated  in  Fig.  3-2.  The 
latter,  W (l-*i),  favors  higher  efficiency  and  weight  while 
the  former  ^ favors  lower  efficiency  and  weight.  The 
difference  arises  from  the  fact  that  as  (1  -A  ) goes 
from  0 to  1;  ri  -goes  from  1 to  oo  . As  the  system  analysis  shows 
p to  be  the  optimization  function  for  minimizing  total 
system  weight,  then  a sub  system  optimized  using  the 
objective  function  C (1  -tv)  will  not  yield  the  desired 
results  but  rather  one  that  favors  higher  efficiencies 
at  the  expence  of  higher  total  system  weight. 

3.1.4  Sample  Run  and  Flow  Chart 

Listing  1 presents  a sample  run.  A user  flow  chart 
is  presented  in  Figure  3-3. 
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Figure  3-3  (continued) 
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Figure  3-3  (continued) 
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Listing  1 
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10. 000 
.1151  E-08 
l . 3262 
* . 1 129 
2. 0000 
? . 00001 


1 05 
1 © i > 
1 8 ? 
118 
1 1 2 
1 1 5 
115 
12(1 

124 

125 
I 26 
127 

154 

1 55 

156 


169 

170 
174 

1 76 
1 8 1 
162 
163 

2 0 o 


234 
tl'  ‘ 5 
245 
2 4 6 

248 
2-1  *:« 
2 8 i..i 

2 8 u 


01  • .. 


22 . 000 
22. 1 8? 8 
2 . GO  GO 
5 8 ii  u 00 
• . .1  000 

1 . 0000 

^ 250OGE 0 . 

. 50000 
0 . 

4<  2-( < 12 

0 * 

0 „ 

2 . 00 00 
0 „ 

1 . 0000 
. € 2*500 
1 . 75 lie 
j 25 . 1 10 
1 7 5 . O 0 
1 . O608 
, 50000 
H 25000 
0 . 

0 . 

. 5 0 00  0 l.  i ' . 
i Oo ooo 
500 00 „ 

, > Ou  i i : • i h • ' • i.O 
. 50000 E U j 

1 i i ■ . . . : . 0 

97 . 000 
97 „ 000 

1 26. 00 

0 . 

0 . 0000 

52 7 . 00 

„ 5u0OnL  0 
: 8 6 „ 00 

0 . 

1 1 0 . 00 


Listing  1 (continued) 

I ' . 1 • ■ • i . i 1 • t • 

E NTER  0 * 0 > 0 TO  END  OPTIHIZl  RAMI  ifl  I INPUT 


0« 

P# 

IN I T I H ■_  V ft L U c. 

1 N L R t r*i  £ 1-1  1 

FINAL  VAtUE 

NUMBER 

i i •! » 

\ 

4b 

1 . 0000 

1 . 0000 

3 . 8000 

0 

' ,1' 

• 0 

48 

2 . 0000 

2 . o 0 0 0 

22 . 000 

0 

I Ci 

: 0 - 1 

22 1 U 

o 

1 09 

16. 000 

1 . 6000 

22 - 00u 

0 

i : o - 

'"*4  i ! ■> 

. 10 

24 , 008 

1 . 0888 

3 „ 060 

0 

U ? 0 - 

ij  Cii'i* 

cr 

0 

0. 

0 . 

0. 

0 

I NF'b  I UP  i I t'i  1 t h r . L'i  I r IJ'NU  ! lUli  DAI  h ■ 

PARAMETER  #•  MIN. MAX. MULT.# 

EM  i Er’  UiU  I u tun  INPUT  uF  UPTIMiZATI l-N  FUN-  T i UU  DHiri! 
:'01  1 


OBJECTIVE  FUNCTION  PARAMETER  NUMBER  = 98  • MULTIPLYING  FACTOR  = 1 . 88 


OBJECTIVE  FUNCTION  PARAMETER  NUMBER  = 0 , MULT : PI.  TING  FACTOR  - 0. 


INPUT  VARIABLE  INDEPENDENT  11  hRhnETEF  RANGE  DhTh  - 

PAR. # j INITIAL  VALUE?  INCREMENT > FINAL  VALUE » #INCREMENTS 

ENTER  0-01 0 >U'0  TO  END  PARAMETER  RANGE  DhI A INPUT  ! 


Ffc 

P# 

INITIAL  VALUE 

INCREMENT 

FINAL  VALUE 

NUMBER  OF 

1 0 ..  r 

t - 6 , . 

1 0 Q . 0 . 0 . 0 

INPUT  ( ONSTRAINI  D B MJNDEB)  DEf  El  Dl  PI  AMI  - R DATA  : 
PARAMETER  It  * BOUND  VALUE  - LIMIT  TYPE 
ENTER  O»0*0  TO  END  CONSTRAINT  VALUE  DATA  INPUT  1 

C#  P#  BOUHDRY  VALUE  LIMIT  TYPE 


1 43  1.7500  -1 

u . ■ 1 

j 0 0 . 0 

1 Ill'll  | PARAME  ILP  LISTING  DATA  : 

Pi  ip.  ills  I EP  k ; ENTER  0 C2ERC  • VO  END  LISIIHG  DA  ill  ilii'UI  ' 


Listing  1 (continued) 

•Mil,.  Dc.'-ICH  i ■ 

Pfti.  E 1 ijp  ' t MiL  ut  : GH  x 
UOtlSTRA  1 hi  D OPTIMAL  DES  Lull 


. . ?••.  -1.:. 


o.SiS  seconds 


11  205  7u 


VI I . 0 rpu  E Ml  ■ ■■  F'HASI 
PRIMARY  LI  HE  INI  1 


L'KI'UL  i '•  l rr‘ 
V’  l\  ^ hL  - Ii1. 


I V • ,:'-'EP hGt. 
h r_. ‘ i 7 ! V i 


r ND  ! 


; 1 r,  I Pll.. »•)  Oh  L y. I i I c:  1 Hf*  V v._‘JuE.N'  i 

HINDI  IG  ( 01 1 E Ui  1 OS  ENSI7Y 
t;.'.!  if*!  , WlNDINu  SUN,  - I UK  .2  ..Ui  I , 
l i ,.l  H-  ■ , r ' 

' i"ii  , 2iJT  Jth-:-  • I 'i 

.iinPY  ht'l'ihu  - 0.  ■ uC.:.  ING  r ATE 

OhjARY  Hi  r ill  i NG  COriL>:JC‘i  OR  CouLIHu  RATE 
• ' 'FT.  ION  KhTL 

, i iF't  LEG  h l i.'TH 
■Ml  LEG  lEPTH 

1 1 PKI  I',,  fRY  i IJP  I IS  r'l  R II  i'ER 

■ i ! 1 f fii  r I !iu  " ACT  OS' 

i i'lAPi  PH'.  :HG  tOUll  ' 1 ■ - 1 
, ..LI  HUM  OF'EPhTu  I'iiij  F:  i i ''LL  S' 
i'Li  E I USUI  nTiOIT  i.*Er  l L , V;  ■ - .6jj  ■ 

fit rORNi  - ' IDfiP.Y  LAVES 

I.II.Hl.  UNDUE  rOr  TEWS'hF'PT  I ip!t 
' 1 ! im*_  ' i.i  i.ii.ii  I UP  TEMPI  i HTL'PL 
. 'Fii  1UKV  I i •HVEPCEHCl.  i.  Pi TE?  iON 
Oil  r;,iF  hit joi ’ CYi.LE 

, , | i.. i . i ■ ■"  ' HI  i • 

ii  , . , JNMEI  . Qi  ■ I FlCT  01  ■ - i iL-  > 

,'L;P  "iilT  FPFAUEHCY  DROOP 

■ WAVE  I ORM  iv  FOR 

i,.'.  i ..  ,L.E;  Isf  ilViNG  .'  Li-  i*stC.  illTE  F'FsiCES 

1 1 ; l . | , , NG  FACTORE’  • 1 " I 

. , 1 0 INVEI  ii'’  ITTERATI 011$ 

. 1 ; 1 1 1 CM  M EH:  2'IENCV 
KH  : 1 1 ui'!  ' 1 ,yi  : EMPtRA  lORF. 
i.iiiii'iP:  i-’i  iWE  S’  • EsP  1 i'GE  KiiCT  iF  iEP 
or  ,.  .Yuh  hi  1GHT 

■ . ,-:i  L r.i'-.  i oLI  Fit  i Eh  F 

I 1 . r ' | ; 0':.  bl'i  ELPOMEH  i 

I I ri  1 os  ■ fplouenc v e::.»L'  1 


■ 'JrO  L OS  ■ 
1 Uhl  1-0''' 

1 I h i | 0G  -. 


4" 

f Kl..r.r' 

■ 

0 . AHt  i l 

•' ' ui 

4 1 

' LCOlL 

TP" 

j i ftriE 

if*,. 

41 

n!  Ium  4. 

HI 

It  TH 

a ; 

ihii  ii 

hi; 

I Gif 

>i  i 

!» 'Ll  P 

or 

PR  3 

MARY  TUI- IG.  1 

a* 

f 'Ji'itjl  1 

1 lr 

stc 

t •tlI<hr,V 

,TF" 

4 

i'UHUlc 

UP 

f-r'  ] 

I1. HP  , t li 

a , 

■ li.'MGLr 

rth 

• 

• 'NIlrtR'i 

,1  'Ll-,  !■'  OF  .1  ' IP  1 H 
1 1 1,  IPt 2.7*  LM  till 

,11  II  . I .1 1 . • I Krill'.  : ' 11  1 


UlYES'  • Ilf 


°t.  COIL 
ft ' ' 1 hVE  !• 


1 00  00  . 

VOLTS'  n'v  i, . 

a-.'  ■ . 06 

VoL  ! $ * rT.  • 

■ . !l  ... 

Ollli  INCHI  • ' 

. 723U2E-0G 

OhM  IUlI-.ES'HV 

90000. 

LINES,  i*  . IN.  > 

50OO . r 

HEF'TL'M::. 

I'.'.iL  1 1 

POUNDS  • L‘ 

- 

. 32400 

PO  1 MI* E CUE  . ■ • 

. '4  000 

PO'  11 1,  ■ ■ S - 

. OUSoOE-Ol 

FOUND'  lit..!'- 

f.7 . 086 

Url  • f ■ : • ' * 6'  1- 

7. 

HrlT  1 'i'I-'LHRE’ 

SE. 30t 

Li'lTT  ■ FOUND 

1 5000 

I Hi  HE  . 

'.  7 SO 00 

1 lift  IE  • 

!7.2i 1 

Ti  IS  1 IS  L-  its.  P 

. 0000 

UNITS  OMITS 

43 

. 140721  -0. 

iw:nt  •: 

1 0060 . 

SIEPT7  • Hu  .. 

sr 

)l  ■ ; HI 

F'lil  .t.  T COO  H 

...  03 

INCHES 

75 . you 

DEG. F. 

: 27 . 00 

DEG. F. 

. i iHjijOE-  02 

.UNIT 

120.00 

■■  ■ .12': 

fit 

G. 

0= i -ROh  _ E 

K 

0. 

tz 

. 5o00i  1 

1 . 0000 

AMPS  fir.i 

2. 0000 

IHTFI-LI  HVIilG 

. ; • 11 4.  02 

IN  < IN -HAT i S 

7 . iioori 

. ?fHi00 

UNIT 

125.  i"’ 

i GF’l 

1 600 U • 

WATTS  • in 

1 . 70  K- 

POUND  ■ 

. 115 1 7E -Gfc 

1 . ‘..if . 

1.  122 

• V5.-4E-U1 

1 NCt 

. hUOOVL  02 

I IN  HI 

ft  >C  . 1 

IIIUHF. ' 

1 .0151 

lN..i:i 

i 7. 261 

TUS'N':  toil.. PS: 

r.-  . • 1 

TURN'.  uOlL'  T. 

1 . UM0U 

LHi’F  F"‘  COil«P 

• 

8".  076 

TUSNS.  LAYIS.S 

•I.I  t I II  10 

1 ii  ,F  s • iu  1 1 « • 

• 7 »6  • 

'UNI  1 

s -1554 

PCI  Ii  ID1 

* 


Listing  1 (continued) 

PAGE  2 OPT  I I’lftL  IiE.--iGH  i 9.515  SECONSS  il  05-76 

CONSTRAINED  OPTIMAL  DESiGN 


c.i 


Ft CULM  I ION 

input  ..ihe  current 
RtC  y ; : . EX'  C HIT  . T CUR,  , 

TOTAL  INPUT  POWER 

Pf  I MARY  WINDING  CURF'I  hT 

r . r'.  ••  i UR/t ,nT  BEUS  . . i > Pi;  . 

>-  J .‘HP  i iivRt  I I I Oil  PP  I < . i'i  - 

R V 1 Q j ID  L - i C B C R 1 ■ E(  F I O N I )L  ( IF  6 H 

! . )U  ■ ? OS 9 -SECT  I ON Ptl  i IREl  t - 

. iGNi.  I : C INBUC  TION  III  ORE 
I 1 r Pill  * T RES  I STANCE 
! GTlii  PR  tip  ' I O- $|  ' . 'I,.',  ■ i I AC  I HE 

• I nil  INC  'J6P  s'  H PRI . +StC.  - ,'ftC  i , G > 

PRI NAP  V i_ Ei IGTH  PER  Cull  ■-  i iLHIl  • 

' P t UHCaP V LEIIGIH  PER  '.;u I L HERN 
I PII'L.PY  RESISTANCE  PEP  60 II. 

S E C 0 NBARY  RE S i S Tfil-JCE  PE R 6 OIL  * 

1)  I S6  i PAT  i ON  PEP  C 0 1 L 1 P'K  l . . PEC . 1 

1 U ’ liL  W . r E VULUME  1 PP  . . r Pm  I . 1 

'.  I 01  NIP.:  WEIGHT  ■ PP i . + SLC . 1 
CORE  VOLUME  , PHYSICAL 
ROPE  WEIGHT  - TOTAL 

■ OKI  p . IPf  i , , INT+I  YSTI  . 1 US  - 

1 1 , 1 . : siPHTi  ccoRi  n ■ 

PRIMARY  COHBUC TOP  VOLUME 
'•ti. unour'i  wipe  volume 
T AUK  INTERNAL  VOLUME 
‘ OQLANT  VOkJJhE 
, UULAMT  WEJuHT 
' RlMAFT  .IMF  PER  TURN 

■ I- L OiiBFiP  r.l'iE  PEP  TURN 

Oil1  I WIKI  0 Pi  i - ’ ,L'  , ll 1 in. . r ,PL  ii 

ll'P-U,  ATE. 'i  MIRE  CkOSULi  T : mNAl  o -i  .u  < 

! . hi  Hi. K IMi.'ii  r.iN'.  E REt  El-V'tli  TO  '.Ei  > •!  ILo  iO'V 
l i oi  i-'.t  PortC inMCt  RK.Oi'T.'  i’U  ■ll.ijI  I To  i i o' 
HOSMucD  IP  r t ML  i r HNP't  UPMI  P DIP  l LA  Err  10 
ABSOLUTE  MINIMUM  SPACING 
T l P Ufi  i PL  Hi  ril.L  L 


Pi  i "till  i Mol.  BAN*  I 
EH  IL  ILTICV  rTP  Puljilli  I'Ll''  1 ! . . i.uiriT  i 
1 Oi'  t Li  G 1 P'P p.  - p.l  L f l ,,'iiHl  I 1 , ' i 
' »..(•  o i EG  MSPLC  r RAT  10  'hi-  f’ ! H W tint 
C 1 1'  ,E  IN  , JPNaL  HEIGH  i' 

*-  ' I.  TIEPTH 
i hSE  width 

' ■LCONIiHP'i  L i I it-  OoRpi  .H 

• i'  Oo  DA;  i'  W .Mil  i No  CUR  P:FN  i 

• i- 1 UNBAR  i.  . I, 

1 • i V I NG  NHVO.i  O' Pit  .•  AG  - I ■ 

I ' ill  '1  CullUHpV  1.1  III  'if  i i 1 ll  t I II 


OU  i 


, P’GSEOE-Ol 

UNIT 

40. 60o 

RMF'c  1 

i . 0000 

AMI  i R!  1 V 

l 1 *“'  . . 

VOLT  - riM'bT'C'P- 

4 0 . 1 l r . 

I li'H  'L  PC  P ' Kl  ' • j P 

65  i .V  . 

AMPERE:  IN 

5 i i 

i 1 il  P E P 1 ‘ ! - 111 

. iS:  . i E ©3 

SQUARE  - 1 ; 1CHES 

. . ■ :4E  04 

SQUARE  INCHES 

0 j » 100  , 

LINES 

. i i 60S 

OHMS.  '01  111 

, 0 ,"0Ci9 

INCHES 

. 1 151 

I ML  HE  P 

ii  i . 1 1 

i Ml  1 HP  C ol l - 1 ' 

4 o05 . 6 

111'  1 i U i L > S 

. 1 

OHMS/ COIL • OP  I 

170. 1 4 

OHMS.  ■:  1.' : L ' SEC 

' • 

WAITS  COIL 

.1  .... 

CUE!  1 1.  I til.  he 

. ! 3£  ,,  | 

i MBS 

. . Oo , 

' BII  .nil 

J . (74 

Pi  LINDS 

0 3 „ 1 o 

hi  t , : 5 - •.  1 •n’C  1 0 

. . , 

NOTTS' '!  RAMP-  T 

. I 699E-01 

CUBIC  - inch;  p< 

. Pi  01 

liJBIC-lllCHEC- 

OP. 095 

i. Ufc:  1 C INC lil'P < 

8. 0507 

CUBII  INCHI 

er  r, 

. Ji_  . _c 

i . . DS»  I 001  AN 

16.5 1.  J' 

VOLT T URN- PR 

1 5 . 9 o-l 

V n,  V-  . TURIN  SE 

0. 

O'UHPL  I HGHES 

0 . 

. UHF . I Nl  iU-I  S 

.922;  L -02 

Hi  NR  I E S 

57fi.  48 

OHMS 

. ] 4 6*  c*  c' 

• UNIT 

0. 

IID.  IE' 

. , r 304 

OHMS  OHM 

. 

OHMS.  OHM 

. , ' 1 

1 .. CP  't  11.  IJlilA 

1 . 1 250 

SOUHPL- INCHES 

„ oiniQfi 

INC  HE'.  INCH 

2.5151 

; in:  i ie  ■ 

2.401  -3 

INCHtP. 

m , o ‘ - m : 

1 1 IT  HE' 

I . UOUA 

ft  iPCPLP ..PM''  1. 

i . U0UU 

HO!  i 1 .1  PP  Hi 

; i 10  m. 

VOiVi.'l  0 ..EC 

. 4 • MU  4 1 1 F • M / 

11 

1 ! '.il  ' , ' . 

VOI  1 I T S ; t 

. 4 • 


I 


i n ,i:  3 OPTIMfiL.  D£':  " ■ 

Coil  FRflll 1ED  OPT  : MflL  Dt  3 1 


Listing  1 (continued) 

OPTIMAL  HE':  TON  i 3.  513 


SEi  UN  1. 3 


■.F  v JtUr-'cE  PKiHHP.  <.  UtsLi... ■ -■ : 'L3  i . H . I 

r"  BE Ci  EE  SI  l ON  BAR'  Oi  I Of  RES*  1 IV 
Of  il'Ul  . E IPcR  ■ f (VI  ?f  ■ 

,3.  '',jHK  ’ ChNThjC  f Ok  ,1  - ' f li  Vr.KMuiS 

R RY  - - 3!  . rOR  El'IPEf  R£  PER  IEN 
£ Curl.  1 EflP't  RJ  ?g  CGEFF  ( RES  S 1 I V 1 I ¥ 
MfiRY  1C-  ■ 10  RE  5TR.NCI 

; -3.  Oi  Di  TC  0 - I R SI 8 i fil  I! !■  1 i 1 , 1 ■ 

' • n. -I'  i 4 1 PE  GAUGE  3’.  L I I H i'JIkl 
Cttui  L-rtR':'  1 1 t:Z  RhOGE  'iO.ul'i.3  3.1  RE 
iUPiBER  Of  IECTIF I EF  B #!•(  E ! I PEE  STfiGR 
/tC'i  ; !;  IER  BRERKIiOWN  c.rlCTuR 
; cl.  ; .3-  ill-  PE  hi:  PEVEPGc  cV3 
.1  1.  , ; fir ' ■ EH i r'E K W .1.  . i 3 1 ! • 'c 
I 1 'i'l'i  c-' 3 1.  fir  DROP  PER  l it1  1 lVIc.1-'  I t ' . ■ 
c i 0 1 i u H i ' V H 

'••PIWhRV  height 
' . : 1 UHJfiP  ■'  HE  I GH  I 
0mPE  TOO!  DIG  RElTE'GTLi  ijH  6 
ERI i 3 I nC  1 1 1 R I HSUlmI 
f i . r 1 1 HUH  RlILFl  I I YE  3>r'H‘.  ■ : NG  .'  rlG  1 He' 

(CLOSURE  HEfiT I N Rfi  E 
. tCOilBOR'i  EH  I PER-  LAYER 
I i N I i'iUl'l  CGGi.  I isC.  GPoCt 
HIT'  i MUM  COOL/ IIG  GRACE 

R IriUl'i  IPERfiTIN  i . 

?fiT  10  C 0 L fit  i 1 VI  IP  OR  1.  I 1 - ‘1  . 

I MhRV-  * O-SET  OHBfiRV  HmF'  1 1 On  1 HI  GPri'  1 m 
EMF  SEC.L-i  4 : Hi  PR 

,1  , Gj  I ■ JIBE  I l¥EF  I ) 1 RE  . 1 1 Li 
Ppfil  . 3 (R  . MMI  If  LfiVER-1  PI  00  53  Jl 

mTh  i MG  • HF'R , VOTER  Lit  i i*  - -El  0.  1.  U'cH,-  i-, 

HI  I nG<  HUH  I .Util  ill.  • 31  1 '.'tie  1 .'  in  HI. . . ' i'l'.'E 

Pi  1 i , .3 1 ,i'  1 1 1 . iii'  TER'  rO  OTEF  fiVI  F 

;ERT  .3  HL  f OP  ' 1 1 . iil.fr’ . ■ I.'  i . Ir'l' 

1 I .1.. Vi  I'l  . til  ■ r I I.II'I  ' i r 1 ,ii  i.<HH  I Hi  ■ 'HI 

, fie  1 hg«  i' IP  ■ . n . 1 ; i . , . ON  1 I ■ 11  IhfiPY 

: IHG» VEPTICfil  FOP  PR3  ■ 1 'i  1 

•Pfifii  INC  > Vtfi  1 1 I T OH  I-  1 1 1 . 1 ) 0 1 

>•  tfitipy  . 1 1 :i ; ; , 1 . Eme 

1 , ’i.  ; 11 . . nO I*  I 301 1 T 1 1 L . Hi  1 l l.  .r  1 ■ * ’ P I I'll  H ' I H 

: RV  ( OHDI 1 1 I OR  ' 0 1 I 

1 in  ,l'  H . rot  ltil.il.  I'M,:  u.r.s  RfiTl.  PL’R  r .3 1 L 

r>  'T;  I .ink  f Mr  t cl'ltil"  i 
t*c  I *•  • hi*'  Ill  111  ' JW  ht.  11.3  - 1’-' 

1 r -Hh'  * 11  1 1 ■ U i‘l  HL  1 l,‘  1 ul*  ‘'I  I* 

ji  Li.  ' • 1'.,  ':.L  1 ONj.ii  iRV  ..  'I.'l  " 

: ' .1  ll'IBEP  1 LtiVi  !•'  1 • 1 1 1 '■  ' 1'  !■  . 

I . '-l!-'  i HijH  MO.  tip  Lh'  .'EHi.  rtr  *'•••  'r 


, SO44Ufc.-0b 

OHM-  i . 0 i t mi . 

= 

, ? 3 4 4 ME  -06 

OHM -IHCHI 

1 0M . UO 

IltGKl  1 r . « 1 . 

r: 

100. O0 

DEGREE  r.'H.' 

- 

. 3 f.72i  IE  -0E 

00(4:1  G-'tE- 

=. 

..  . ' 

I. 

= 

i . 0080 

OH  1 1 '.'Hi  1 

ss 

1 . 0032 

OHMS  OHM 

5= 

8 • . 0.00 

S T PH  OI10 

= 

0 . 0 0 0 

S I'HhH  n: 

= 

4:3, 521 

RECTI FI  • 

= 

' . 00  00 

ViJLJ  3 v.- 

s$ 

5 D 0 . Q 0 

vol  rs,i  1 

=: 

1 3 100. 

v 0 n Hi  n 11  ; 

1 . 1 GH  11 

’•.  lil.  i ■'  . . i- 1 ILIlT 

1 1100. 

VO  ••  ‘ 

,51249 

1 til-HT 3 > 0 R 1 . 

m y y h , 

IHCHLS  E'  . 

=r 

-:l  1 7 bo 

WATTS  50.  H. 

i „ 0000 

VOL  1 S VOL 

0 . 

I II CHE S INCH 

fl  . 

WfiTTS  30'.  1 

1437 . 8 

VOL ! ' LAYER*  1 

•ss 

. £?8OOE-0i 

U It*  11.  • 

X 

. 80000 

i iliC  HE  1 

ss 

0. 

PHD  i fit  i3 

X 

. 4 ’.Vu HE- 02 

(j.-ilj,  Im.  4 

3= 

. 15665 

INCHES 

0 . 

VOLTS 

. ' S4SS 

INC  Ilf ' 

. 1 654  0 

INCHES 

. 1 j 8 9 4 

i he  he"- 

a 

. lH5H4L-tH 

111!  HI  ■ 

s: 

. 12286 

i 111  hi  . 

s 

. 2550 i 

INCHES 

a 

. H6HHuE-01 

I I IC  HE  :• 

s 

. 12144 

1 hi.  H ll 

58 

. 24  0 1 

in.;  ut 

a 

.1. 1 '.*  4 .L  i.i  i i m:  Hi 

inj 

VOL  1 .'IOC  «PR  1 

.HH3iH:E-m1 

INCHL3 

219. 38 

WAT  To 

1 70.  14 

Ni  l r 1 5 

6 O 3 • 9 4 

U H 14 1 • 3LI  . 

. . * 

I liC ML'  -l  l-  1 . 

, e.t./:*S8 

INC  HE  • 1 . 

2 . » 1 U i '11 1 

GROUPS  EC. 

4 . t'loeo 

0 0 HI-  i.RU'-il'  r 

5 . uOOi » 

LfiYEF  IROUPil 

. ' „ nm.M  1 

i.ROUf  :■  ' i-i  . *r 

Listing  1 (continued) 

P H ‘j C_  4 OPTIMAL  DE-mGN  i »,.5l5  SECONDS  11/05  -•  70 

'•  ‘N  ■ - PH  : HE D Jr  T L Mi L JL  7 < i 


...  GROUP!  SEC,  • l 

!LhH  COIL  PERIMETER  ^ 6. 61 13  I NO  HE': 

V£  €01  ; ■ SQUBI  HI  . 

1 -r  h E~  ro-  ‘ HSE  HERD  oPRl  t*  Tu  hL  = . 9000  j 1 Nf  Ht5 

• f:?-  •.tL»nOf';._  iMTiiL  ' • PuR  7 STf  U»  • *r  UPE  - 0.  I NChE  :• 

I • U i KI i i\i»  C.’Eh:  IiL4‘.  i F Fit  , - jP  = 2,i  UU  VULTS L.T 

1ARY  ' „ • • • GRQUI 

. . ' ' • ' i • , ■ 

. . M i * _* 1 1 1 ! Or  LH  i « Pfc  P ' P * i • ‘ah’ • V'  " 2.0000  LhVER  OOPiM 

o.  ri  I 1 1 1'iUt'i  LH  EP  PP  i PfnF’V  LM'IO  = 1 . QkJkfQ  GROUPS-'  l b:  i ■ Mi 

• . iJc.i  ' MhLirlU  i LflVCF  !-V i MARY  GROUPS  - U.  GROUPS  rfcl . .«  M 


. ■ n .l  t:  . L . ' U“ • OL  i i* 
u ir  L Hire.  1'IRf‘jETER 

■ < ' 1 • H - i.‘  5 . KE  J I RHE  . ER 

l"  . 7 V F 1 1 »_  . UR 
. I . 'F  iGHT 

. • *5  ..  . r , r I . HIE.  , 

■ ’ • r lUN  F* ri ; H DEF‘  i H 

UHiiETER 

'.'f  f,H*  . li'j  l-'  i Li  7 nL l'EP  i ri.' 

. UMfcLR  i.  . ROMS  DEEP 
' <i ?!-.’•  Ur  r'Oi  ’ 7 H • EE 

• i i . n ’*•  5 FLUil  RhTF 

rlWl-.ij,  i . !•*  ' PMF'ERRT1 1 -’t 

i F 'ii.it  •.  •-•'.4.  I IN  ; » EUr’LKrl  I '.'Pi:. 

• ! r Hr'1  I 7 P tt  N r.  F F F COeFF  I •.  i 1 i 
Hi  i 7 .iBUCEp  SURFRCL  HP  2 H 70  Ml!' 

riH  i.1  . h f i t I (' 

ii£ r r»» t ; , i.ajt  mkhli-  t hEFO  i ,jU!, 

OLU ''•!*.  7 . ’Rm  .T  URl’L.t*  + HEh  • *■  ! , 

*L*  i HL  il E '.CFlrih  iLt  LL  . **•  *■  1 . 

'.'"IHHT  "T’.'F  T'UPL  DROP  <i.r.  £/!...  r,;IM7 

• ' i *L  I ..  I .Oh  r»  H , t-  ij.  f . E l i . , 

F 1 * I .1  Our* . Kt  • 

Ml*  ' * *-  -•£*-  . *.  «JhD . r’ES  . a v 1 . 

■'  ij-  F‘P)  J uNIi.R'i,  i f .Mp.  - •»  “ . i JL-; 

1 « • l nL  ' OH!1 , f l.  . I . i r . nl’H  . * • • ■ 1 

. r ; *-mi m .i-  f ':«*  rfb*.  r.  • . ; t j-’  . ■ > 

• • ihLr  '•/.IPUP  • -ul  •*  ..  .i  •*  . • iOil  •.  i.*i  '1 
u:  1 4'iLP  H p .L  r OHER 
» *.  l t •EoF'ERHruRE 
« *■  7 • HHMI.LR  Hj.  i GH 

■ I '•  :JPl  JF’OP 
f i P VELQF IT  7 

• l4t.r‘  i in"! 

. ..JU  W;  1 IjH T 

i . .J.'  (h.  t<  if.  ’ (*f  . f t 

t-’Hi . i !t>'  .]  u-,  |.r  r.*  i • 


c = 

0. 

• = 

0 , 

I NCIitS 

Ei . 

i ; i.‘  he 

4 . 0000 

'i.ECGii  - i ■ 'icCON 

0. 

POUHL'S 

0, 

lll-.  l-: 

tl . 

I HCHES 

. 6250O 

1 HL  ME  .• 

1 . 7500 

1 HCHES 

0. 

ROUS 

O. 

ROW 

= 

0 . 

LB'.:  v Li  ' . i M 

J 25 « 00 

il  Ij  . [ r : r . Ii  ,i  . i.i 

- 

1 75.00 

IltliKE i L I iKt.  I t 

= 

O . 

JOUL  ' I.  . I 

0 . 

snuopl:  Ilf:  111  .• 

o . 

i m: ; 

= 

0, 

INC  HE;. 

0 . 

CUBIC  1 i 'lit, u 

1 . 0000 

CUBIC  i . : ■■ 

= 

, 50000 

. 25O00 

. . • • ■ \6 

OHM  ....  II  • 

, , ...  | * ; - 

ii.  t I NCI 

.... 

. 

. a 

, 2U138.E-02 

DEC. . . 0 

0 = 

0. 

SC-'I  l| . , III  Hi  v. 

ic 

U . 

jijUL  '..Li.  ..I*.  1 

0 , 

Nil  f i ■ 

IC  .•.!•'(  ' ■ . . i TN 

I Ni.  lie 

0 . 

III'  lit  c.  1 . nil 

y, 

t mi  Me  j. 

0 , 

i Ui  |I  1 H . ill 

0 * 

POOL  i 

1 : 'LH  U' 

i'li.  '1  .i 

U • 

PCliIM 

J — i 


Listing  1 (continued) 


F REE  • UPTiWRL  tiE'i  IGN  i 

C 1 L' 

'Mi  COHIi'.- 

1 1/85/76 

-■!  r.-Ti  'ii ; NEIl  UET  liltiL  Dl'cUUh 

v u.  Li*.’ ! i N'"*  T , -Ji  i c1»lL  jUKFmL  L P.REft 
. _ .-f i . ni  _•*. -k'Fftb : HR£H 
1 E l Hi'i  i 1 ■: . uh  L I*IjE  ._:uGL  Pn  1 1 

Fft'Ci  > ■ ■ i 

■ _ . , i'i  .ii.  '.'.‘I  n I"-.L  FftCir.  L-UUL  KH  rL  Phi  ’ 

OMBi  IN  CKKESS 

• vr  n.’L  j.  rib  Us  .rlT  I *.Jl  -i  I H ».  Pi-  . itob 

• < . U.  ■ • ■ I L. 

• • • ■ GE  E R - 1 IG 
. . r.;  , R ft  t.  c K DIE  L E C T X « 3 ; r ! EH  G T H 

ft  . COO  ' ; 

■ f f ' CON  FftCI  . Of  •;  - . . fOR. 

■ !•;  Jr* . L U.  U-  Ffi  .-c  I.  UUl  hr  rk-r  VLR  I I H I L i:. 

•J ! L Fuki’i  INNER  FfiL-E  fN  iPKiLioG 
1 h.L  FOrtN  OUTER  FftCE  TH I ; t!E:,P 


' * 1 

I.  f*  •<  L ' M JER  j”  I :bh  ' i ; . L : if  • :• 


i ~'-r 

_ ft 

r IN 

r.  ij . 

BARRIE R fHICI 

HESS 

r ” *_  P ' 

I LH 

,'ER  IN 

;.ij . 

EhPR  i Li\  rUl  .-l’ 

HESS 

; R PF 

: Lft 

r'£R  I N 

EJ . 

, ER  , i i i 

0 

Y £R- 

F.NCL0 

INSL 

OPlRP  iEi'1  i H I » 

NESS 

( 1 - ;■!  , 

i L n 

V INSU 

. 6RF 

• . Lft  ’’H  [Ef  MLS 

5 

■:il  hl 

CuO:. 

SPACE 

TH  l L f F i.'K  ij  El  Eu 

: i-AVEP 

r o > hL 

COOt. 

Cl  '.-ICE 

THi 

. 1.  FOR  ONE  PF: 

i LAVER 

j.  fiOL  •£■ 

POEE 

'1  i.'O  I RE  h 

hT  TOP  O'-  PR  I 

LAVER 

\\  . 

PftCE 

P'E'.iUl 

' C I1 

AT  TOI  OF  o - 

AYE . 

• Ui'L 

1 1 . 

PEOUI 

RED 

Ml  BO i OF  PR I 

LAYER 

COCL  . I i 0 I fiV-cF 

I ' -E  ETiifi  UP  ft  NT*  LOi,  • .'j  IL  Fu?-f‘J  r dCE 

1 1 FI4N  U BN  1 FORI  FftCI 

, . , . i . , i . . BftR  2 

i . N L i 1 0 . L ;jh  ! ' . PEL  j F I »_•  Hr  hi 

. 

p ORE  OPEC l F [ C HEftT 

. NB.i  OHB.  H R HG  PEI! 

I . • . i . • . OHD*  IN  i ii  . , , . ■ . 

• .i  'i  i i1  FRftNSIEN  RftTl 

W,  INI).  : I 

• ' .iUE.il.  P Puli'.  OFF  ( ■; . TkkOO  . Un.iL  : Mb  F.rO 

.iPF  OFF  • IMP  TRRM’r-j!  Hi  LEi*l;N*j  F-’ftTL 

. 

H!  EE*  .EuHD.V0L.T0  GQl't  IMG  OORFREE  ;:.R 
• uel  . OHl  lb  ' ULb  : i{ij  EijPFh*  F r;PE 
EHJOf  i YCLE  Ur  F T J lE 

i •i'-LhlU  EE  tV  URL  Dr  • ■ , IH.T  1 ..  MlNiU  NGG 

G viObr,,.  T ; PEG  triC  HERT 
i V h V * »FE  Br-I  JR  I R i 1 • T J .*<  I P lit  J OK 
, ',KiT7 1 ; fiT  IOIihI  ' uNOrh'IT  -* 


= 

e. 

So. IN. 

= 

0. 

SO. 3 h . 

0 . 

WATTS  '.0. 

IN. 

= 

0 . 

WAfT.  SC. 

IN. 

SS 

0. 

! HCHiC 

- 

0 . 

INCHES 

= 

. * UUUO 

1 1 1 L HE  h 

I ii  0000 

I V-  JH. 

X. 

50000 . 

I-:V-  lil. 

5E 

. 280 C ■ 1 - ‘ 

K V.-  1 rl . 

= 

0 . 

INC  HE'-. 

. 6 o E.  Pi  PE"  Pi  3 

INC HE S 

0 . 

INCHES 

= 

. 6 488.  : 01 

INCHES 

5= 

. o4  E I UE*~U  1 

INC HE' 

. i:.H05E-O2 

INCHE 

= 

« t>  IJ  o E 4 1 U i 

INC  n:  ■ 

= 

...Cl.  51  -01 

INCHE': 

, ] Eh  j 

INCHES 

~ 

0. 

IilCHE'- 

a 

. 11353 

INCHES 

. 1120< 

INC  HE'S 

.54084E-02 

i NCHES 

- 

. 51  3t  : 01 

INCHE 

a. 

. . ' 

INC m.  ■ 

. 2!  • 36' 

IHChEs 

= 

. 1 ;. 95E-02 

INCHE 

. 15516E-01 

INCHE ‘ 

a 

. E'E'i'P? 

] NC  HF 

= 

.81511 

INCHE' 

= 

. 1 0429 

In  h.  ■ 

a 

3 E . 000 

JOULi  '•  1 P 

I' LG 

a 

07 . Pi  hi  i 

JOULE'  It: 

1'ISC 

= 

I j •.  60 

JON  US  IB 

I'F  G 

G 6 . 7E0 

Hi  ITT!  O. 

IN 

a 

?.  :0?V 

MAT i .MO. 

11^ 

-1.1911 

NAVI'S  u- ' . 

i)l 

. • . 61  01 

On)  * ' • m.  , 

. ; ' ■ 1 02 

lull  .N. 

IN 

a 

,7  i z 

HAT!  n. 

m 

a 

19E  "1 

cu.m 

IN 

a 

. 4 m.-  .4.1  E - 0P‘ 

Cl.  Ill  S'1. 

IN 

a 

0. 

1 U.  Ill 

IN 

90  Pi . P»h 

SEC  ONES 

a 

. M00 0 

LB',  ' . ill 

a 

5 J . 00 

JOIJL  L . Lt 

DEC 

a 

.500801  -ii 

. . MM 

: in  h1  ■ 

-.c  s 

Listing  1 (continued) 

PAGE  6 OPTIMAL  DESIGN  1 . 

1 3 

StCONDS 

l/@5  ■ ?t 

t.  uN'.  TRAINED  OPTIMAL  DESIGN 

2 j i 

'Ki.77  I'liHDifiL  L-OOLAN  1 VELu.  T |V 

0 . 

INCHES.-  SEC 

■4.  •_*  *_ 

; 

0. 

I NCI  IE  .PC 

1 ‘ 

. . 

= 

0 . 

I NL  pit'  - SEi 

C7  / 

7 l ! i pi  I is  J . COOL  AH  i MAS':-  •*  Rh  i u 

= 

0 . 

NLBS- SL  C.-'IN 

r*  r ; 

L7  - RV  74  i Nil . LUOLfiN  l MAS '3  LOU  E.h  1 t. 

•-= 

0 * 

LBS/:-  F IN 

2-  ■ ■ : 

: is  LUi 

= 

0 . 

LBS.-  . ..  In 

c-  -- 

r r : ;-!r  ' 0 i N.!J . LUULAN  i • Ef'Ir'L  KM  1 URL  K i ' -L 

= 

0 . 

BEG.F 

"i  r •_ 

iR’i  1*4 1 ND - C OOLftNI  ,r  MPERATURE  RISE 

= 

0 . 

DEG.  F 

c; . j 

‘■L'P'b.  LUULhMi  TEMPEPB  1 UKE  F i S c 

= 

0. 

BEG.F 

2*  k 

HCOI  IMG  C 301  ftH r TEMPI 

= 

110. 00 

DEG.  1-' 

2 6 1 

G.1 1 i 1 1_.  L DOLAN  i t EMPERn  : CPE 

- 

0 . 

DEG.  F 

.it 

COOlRMT  THERMAL  CONDUCTIVITY 

= 

. I83O0E-02 

ji’i,  iL.EC  ' L*  •.  - SO 

■ i 

i 1 u L.  A M T V I S C U 3 1 T Y B T Y UP,  IT.  t: 

. 69400E-04 

JLBS.'Il  SEC 

2.64 

l 'Ll  COOLANT  VISCOSITY 

, 54 . i . ; 

JLBS  • IN  c 

PP ; . : ; l ND  -COOL. HEAT  TRANSFER  COEFFICIENT 

=S 

0. 

JOULES  SEC-  SC! 

l 

. M INB  -C(  OL • HEA1  ! F fit i >1  ER  ! Jt  FFIl  - El  1 i 

- 

0. 

JOULE-  SEC  Su. 

26 

■ UiPE-l uOLANl  HER!  TEMNSFEF  LUEFFILIENT 

= 

G . 

JOULE  SEC  SO. 

r. 

NECESSARY  PRIMARY  MINDING  COOLING  RATE 

= 

6G . 780 

WATTS'CU  ill 

■ 

i *r  U r.b  SUK\  >,ELONDHRY  WINBiNC  CUNDUCTuR  UU 

= 

7.307? 

OATTS  Si..  IN 

270 

NECE  iftRY  CORE  COOLING  RATE 

- 

. 55tuj9 

WATTS  : -U . 1 N 

27 1 

1 MARY  WIND. COOLANT  MASS  1 LOW  Rfi  E 

= 

0. 

L P S-  ' • : 1 

' ti  ON  JAR  Y WINDING  COOLANT  MASS  FLOW  RriTE 

= 

0 . 

LBS -SEC 

, UPE  L-UUl.hNI  MASS  flow  pate 

0 . 

LBS-'  SL,. 

2 ? 4 

TOTAL  COOLANT  MASS  FLOW  RATE 

= 

0. 

LBS/SEC 

. 

EFFORT  4 FLOW  - PRODUCT 

= 

y*. 

2 T v 

UP  IllriiYT  COOLANT  FLOW  RENBLDS  NUMBER 

= 

0 . 

L >'  i 

•_.Ei  UNuAPY  COOLANT  FLOW  REi 'ALL'S  NUiTc-EP* 

~ 

J . 

273 

i OPE  LEG  COOLANT  FLOW  PEN  ALT'S  NUMBER: 

= 

0. 

279 

F.jpi.f.D  COOLING  CHANEl.  SPACING  INCREMENT 

= 

. 5U0WWE-U2 

'•■TOP  NORMAL  JOB  EMECUT  i ON  COMPLETED 
1.322  CP  SECONDS  EXECUTION  TIME 

I 


OF  H 12.250  SEC.  4.407  BDJ. 

CPUS  9.258 

C OHNEC  T i ME  o HRS*  34  Hit  . 

11  05--  ?t  LOGGED  OUT  AT  12.  1C.  33. 


i 


3.1.5  Recommendations  for  Improvements  to  Present 
Program 

There  are  two  general  areas  in  which  the  present 
program  can  be  improved.  They  are:  1,  the  transformer 
design  procedure,  and  2,  the  central  and  optimization 
procedure.  Improvements  in  the  transformer  design 
procedure  to  incorporate  appropriate  physical  and 
thermo-electro-magnetic  relationships  to  design  transformers 
having  structures  not  presently  incorporated  in  the  procedure 
(i.e.,  PIE  windings,  etc.)  and  allow  for  a more  general 
selection  of  independent  and  dependent  parameters  by  the 
user  are  recommended.  Also  recommended  are  improvements 
to  the  control  data  I/O,  and  optimization  procedures  to 
increase  user  and  program  efficiency  and  ease  of  utilization 
and  , reduce  execution  times. 

3.2  10  KW  Transformer  Development 

3.2.1  Design  of  Original  10  KW  Unit 

Numerous  parametric  sensitivity  analyses  were  conducted 
to  derive  the  design  of  a 10  KHz,  10  KVA  inverter  transformer, 
with  both  ferrite  and  Orthonol  cores.  Table  I is  a 
comparative  listing  of  the  parameters  for  optimized 
transformers  with  ferrite  and  Orthonol  cores.  Figures  3-4  & 3-5 
are  sketches  of  these  optimized  transformers. 

Significant  comparisons  are:  (1)  Weight  (K50) ; 0.909  lbs 
Orthonol  core  vs  2.47  ferrite  core;  (2)  efficiency  (K49) ; 

93.6%  Orthonol  core  vs  81.0%  ferrite  core;  and  (3)  total 
volume  (K77);  12.92  cu.in.  Orthonol  core  vs  31.24  cu.in. 
ferrite  core.  The  weights  and  volumes  do  not  include  any 
auxiliary  heat  transfer  apparatus.  If  this  were  included, 
the  Orthonol  core  transformer  would  show  a further  advantage 
in  weight,  due  to  its  higher  efficiency. 

It  is  interesting  to  note  that  the  core  dissipation  (K73) 
is  lower  for  the  ferrite  (8.78  watts)  than  for  the  Orthonol 
(23.3  watts),  but  that  this  advantage  is  insignificant 
compared  to  the  copper  losses  (K68)  of  2335  watts  for  the 
ferrite  and  659  watts  for  the  Orthonol. 
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TABLE  1 
10  KHZ 


f ORTHO NAL  FERRITE 


1 

OUTPUT  VOLTAGE- VOUT 

1 0000 

1 0000 

2 

INPUT  VOLTAGE-  VIN 

285 

£1  cr 

PRIMARY  RESISTIVITY  - PP) OHM-IN 

6. 7 8 E~7 

6.  ?&'£-? 

4 

S E C 0 N BAR V R E S ISTI V I T Y - P S - 0 H M - 1 N 

6.78E-7 

6 . 7 y E ~ 7 

er 

F'EH  K F L U X D E N S I T Y > B M A X > L.  I N F S / 8 Q .IN. 

96775 

16130 

6 

FREQUENCY) HZ 

1 0000 

1 0000 

■J7 

PRIMARY  WIRE  DENSITY- LBS/CU.  IN. 

. 324 

. 324 

O 

S E C 0 N D A R Y W IRE  D E N S I T Y , L B S / C U . I H . 

. 324 

. 324 

9 

CORE  DENSITY) LBS/CU.  IN. 

. 296 

. 276 

10 

COOLANT  DENSITY) LBS/CU.  IN. 

. 0505 

. 0505 

1 1 

PR  I MARY  HEATI NG  RATE - HRP - WATTS  -"SQ . I H . 

11.8 

14.5 

12 

S E C 0 N II A R Y HEATING  RATE)  H R 3 ) W A T T S / S 0 . I N . 

2 . 2 

•”j  cr 

13 

C 0 RE  Ii  I S S IPATI 0 N RATE)  P C - W A T T S / L B . 

130 

li 

14 

CORE  LEG  WIDTH- A. IN. 

. 3 1 6 

. 548 

15 

CORE  LEG  DEPTH.- B)  IN. 

. 316 

. 548 

16 

PRIMARY  TURNS/LAYER) NDP 

36. 2 

6 9 . 8 

17 

C 0 U P L I N G C 0 E F F I C I E N T ) C C 

. 999 

. 999 

18 

PRIMARY  SPACE  BETWEEN  TURNS) SIP- IN. 

. 001 

. 001 

19 

P R I M A R Y S p A C E B E T W E E N L A V E R S - S 2 P - I N . 

. 05 

. 052 

28 

P R I M A R Y -TO- C 0 REE N D S P A C E < V E R T > - S 3 P -IN. 

. 065 

.•  168 

21 

S E C 0 N D A R Y S P ACE  BET  W E E N T I.J  R N S - S 1 S ) I N . 

. 001 

. 00  1 

2 2 

SE C 0 N H A R Y S P A C E BET W E E N L A Y E R S ) S 2 S -IN. 

. 05 

. 05 

2 2 

S E C 0 N H A R Y -TO- C 0 RE  E N D S P A C E < V E R T ) - S 3 S - I H . 

. 133 

. 241 

24 

EFFICIEN C Y I N C R < M A X . A L L .ERR 0 R > - D E C I M A 1... 

. 000  1 

. 000  1 

■i  cr 

C 0 0 L A N T HE  A T 0 F V A F'  0 R I Z A T I 0 N - H F G - B T IJ  / L B 

63.  12 

63.  12 

26 

F'  R I M A R Y C 0 N N ECU  0 N F A C T 0 R 

1 

1 

il  (' 

S E C 0 N H A R Y C 0 N N ECTI 0 N FACT 0 R 

1 

1 

2 o 

P R I M A R Y WIRE  T H C K < R ECT. W I R E 0 N L Y > 

29 

S E C 0 N Ii  A R V WIRE  T H C K < R ECT.  W I R E 0 H L Y > 

- 

30 

NUMBER  OF  INTERLEAVINGS 

4 

4 

31 

C 0 0 LING  SPA  C I N G F A C T 0 R -IN.  < WATT  S / S Q .IN. ) 

1 . 53E-3 

1 . 5 3 E - 3 

F'RIMARY-TO-CORE  END  SPACECHQRIZ) ) SPC-  IN. 

. 05 

. 05 

”■  O 

P R I M A R Y -TO-  S E C 0 N D A R Y S P ACE-  S F'  S -IN. 

. 05 

. 05 

34 

SEC-TO-SEC  SPACE (BETWEEN  COILS) - SSS- IN. 

. 05 

. 05 

3S 

OUTPUT  POWER- WATTS 

1 0000 

1 0000 

36 

AIR  GAP  LENGTH- IN. 

. 001 

. 001 

37 

RELATIVE  PERMEABILITY 

15000 

5000 

38 

EQUIVALENT  CORE  LENGTH  - OUTER  LEG- IN. 

- 

- 

39 

EQUIVALENT  CORE  LENGTH- INNER  LEG- IN. 

- 

- 

40 

PRIMARY  WIRE  DIAMETER) DIAP- IN. 

. 032 < 20) 

. 022  (.  20 

41 

SECONDARY  W I RE  D I AMETER - D I AS -IN. 

. 005(36) 

.004 < 38 

42 

WINDOW  WIDTH- E- IN. 

1 . 039 

1 . 032 

43 

WINDOW  HEIGHT  - D- IN. 

1 . 359 

2 . 6 9 2 

44 

NUMBER  OF  PRIMARY  TURNS'.  CO I L. - NP 

-?  •*.  o jZj 

139.64 

45 

NUMBER  OF  SECONDARY  TURNS/CO  I L - NS 

2709. 9 

6049, 4 

46 

NUMBER  OF  PRIMARY  LAYERS/ CO  I L - NF'L 

2 

2 

47 

N U M BER  0 F S E C 0 N D A R Y T U R f I S / L A Y E R - N S Ii 

182.2 

419.9 

48 

HUMBER  OF  SECONDARY  LAYERS/COIL  - NSL 

15 

15 

49 

EFFICIENCY 

93 . 6 

81.0 

L_i_ 
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II 


TABLE  1*C0NT 


I 

50 

TOTAL  WE I GHT - WTOT * LBS . 

orthonal 
. 909 

FERRITE 
2.  47 

51 

REGULATION* REG 

93.  7 

81.0 

cr  ~* 

COOLANT  FLOW  FATE*  LE'  SEC. 

.0102 

. 0352 

C j •-« 

0 U T P U T C IJ R RENT  , 1 0 U T * H M PE 

1 

1 

54 

INPUT  POWER* PIN- WATTS 

10681 . 4 

12342. 3 

cr  it 
. j . j 

INPUT  CURRENT* II N* AMPS 

37. 48 

43.31 

56 

PRIMARY  CURRENT  DENSITY*  JF' * AMPS'  30.  IN. 

46656 

51414 

cr  “? 

j |- 

SECONDARY  CURRENT  DENSITY* JS* AMI  3.  SQ. IN. 

50957 

6 9 9 9 5 

58 

PR  I MARY  W I RE  CSA  * AP  * SQ .IN. 

8. 03E-4 

8 ! 42E-4 

59 

S E C 0 N II A R Y W I R E C S A * A S * S 0 .IN. 

. 196E-4 

. 143E-4 

68 

FLUX  IN  CORE* LINES 

967 4 . 8 

4838. 7 

61 

PRIMARY  DEPTH*  DP* IN. 

.114 

.118 

62 

SECONDARY  DEPTH* DS* IN. 

-7-jcr 
■ i i J 

. 764 

63 

TOTAL  W I HD  I NG  DEPTH  * DSP -IN. 

1 . 039 

1 . 032 

64 

P R I M A R Y LENGTH  C OIL*  LP  * I H . 

153.5 

427.5 

65 

S E C 0 N D A R Y LENGTH / C 0 1 L * L S * I N . 

13800.5 

36188.5 

66 

PR  I MARY  RES  I STANCE  CO I L * RP  * OHMS 

. 13 

. 34 

67 

SECONDARY  RES  I STAHL  E CO  I L » RS * OHMS 

476.8 

1696. 9 

6 8 

It  I S3 1 PAT  I ON - CO  I L * P J * WATTS 

658.  8 

2334 . 8 

69 

TOTAL  W I RE  VOLUME  * VW * CU .IN. 

. 394 

i Oi  ( 

70 

TOTAL  WIRE  WE I GHT * WW * LBS . 

. 128 

.284 

71 

CORE  VOLUME* VC* CU. IN. 

. 606 

2.891 

r <£. 

CORE  WEIGHT* WC* LBS. 

. 179 

. 798 

73 

C 0 RE  D I S S I P A T I ON  * P C * W A T T S 

>•  ~j  3 

8. 78 

74 

TOTAL  D I SS I PAT I ON  * PL  * WATTS 

682.  1 

2343.6 

75 

P R I M A R Y W I R E V 0 L IJ  M E / C OIL*  V W F * C U .IN. 

. 123 

. 360 

7 6 

SECONDARY  W I RE  VOLUME/CO  I L * VWS * CU .IN. 

. 270 

.517 

77 

TANK  VOLUME* VT*OJ. IN. 

12.  92 

31.24 

7 

COOLANT  VOLUME*  VF*  OJ. IN. 

11.92 

27.47 

79 

COOLANT  WE  I GHT -WF* LBS 

. 602 

1 . 39 

80 

PRIMARY  VOLTS/TURN 

3.  94 

2.04 

81 

SECONDARY  VOLTS.  TURN 

3 . 6 9 

1 . 65 

y y 

PRIMARY  SELF  INDUCTANCE* CENTER  LEG*  HENRYS 

13.9E-3 

95.3E-3 

83 

PRIMARY  SELF  INDUCTANCE* OUTER  LEG  * HENRYS 

- 

- 

84 

SEC  SELF  INDUCTANCE* CENTER  LEG- HENRYS 

19.5 

178.8 

85 

SEC  SELF  INDUCTANCE* OUTER  LEG* HENRYS 

- 

- 

86 

I NC REM ENT AL  EFF I C I ENCY - DEC  I MAL 

- 

- 

87 

MINIMUM  SPACING* IN. 

. 05 

. 05 

LEAKAGE  INDUCTANCE  PEL  TO  SECONDARY*  HENRY'S 

. 015E-3 

. 037E-3 

89 

LEAKAGE  INDUCTANCE  REL  TO  PRIMARY, HENRYS 

24. 1 E — 3 

73. 9E-3 

90 

FIGURE  OF  MERIT 

- 

- 

91 

CORE  CSA?  SQ. IN. 

0.  1 

0.  3 

92 

CORE  ASPECT  RATIO=DEPTH  WIDTH 

1 

1 

93 

TANK  WIDTH* IN. 

2.71 

3.16 

94 

TANK  DEPTH* IN. 

2.39 

2.61 
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SCALE  : FL//  L 
SP.A/T.  0.25  LP/ KVA 


10KVA  lOKHz 
INVERTER  TRANSFORMER 
CORE/COIL  ASSY. 
FERRITE  CORE 


It  appears  from  this  data  that  orthonol  cores  show 
significant  advantages  over  ferrite  cores  for  this 
application.  Additional  data  will  be  generated  for  10  KVA, 

5 KHz  transformers,  and  for  200  KVA  transformers  with  heat 
exchangers . 

All  major  parts  and  materials  for  fabrication  of  the 
10  KW  transformer/rectifier  assembly  are  on  hand.  The 
primary  coils  will  be  wound  with  f 20  AWG  heavy  Formvar 
wire.  A special  double  coated  wire  by  Viking  Corp.  which 
has  the  property  of  self  bonding  has  been  chosen  for  the 
secondary  windings.  The  mandrel  parts  have  been  submitted 
to  U.B.  for  fabrication. 

Detailed  specification  of  the  individual  windinas  had 
been  with  held  pending  the  outcome  of  the  review  meeting  at 
AFAPL.  Following  the  meeting,  additional  tests  were  made  of 
the  Orthonol  core  which  showed  that  all  instabilities  and 
distortion  are  eliminated,  if  excitation  is  kept  below  14KG. 
Similar  tests  were  run  on  a toroidal  Supermalloy  core  of 
approximately  the  same  dimensions  as  the  Orthonol  core,  which 
showed  that  the  excitation  must  be  less  than  8KG  for  this 
material.  It  is  expected  that  introduction  of  an  air  gap 
would  force  a further  reduction  in  excitation.  The  maximum 
flux  density  is  shown  plotted  as  a function  of  RMS  volts 
per  turn  for  the  Orthonol  core  in  Fig.  3-6.  The  expected 
induced  voltage  will  be  3.2  volts  per  turn.  This  was 
experimentally  verified  before  proceeding  with  the  winding 
specification. 

3.2.2  New  Optimized  Designs 

Using  the  CDC6600  TD0P2  transformer  design  program, 
a number  of  new  transformers  were  designed.  These  have  a 
shell  construction  using  two  "C"  cores  as  shown  in  Fig.  3-7. 
Materials  considered  included  both  Orthonol  and  ferrite. 
Cooling  was  by  Freon  113  vaporization  in  all  cases.  Both 
special  and  standard  cores  were  considered.  It  was  found 
that  several  standard  cores  could  be  used  without  a sub- 
stantial sacrifice  in  overall  weight.  Some  of  the  more 
interesting  designs  are  listed  in  Table2.  These  are  all 
for  layer  wound  coils.  A program  modification  to  permit 
consideration  of  pie  wound  coils  in  this  program  would  be 
desirable  prior  to  designing  the  200  KW  unit. 
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/ c distortion  threshold 


V/N  - Volts  per  turn 


Maximum  Flux  Density-vs- Volts  Per  Turn  For  Orthonol  Core 

Fig.  3-6 
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TABLE  II 
Design  Data 


Material 

Part  No. 

a 

b 

d 

e 

WT.Lb. 

eff . 

Orthonol 

Special 

.775 

.775 

1.205 

.802 

1.43 

94 

Ferrite 

Special 

1.323 

1.323 

1.892 

1.333 

5.08 

92 

Ferrite 

Special 

1.414 

1.414 

2.912 

0.957 

5.86 

94 

Orthonol 

Special 

0.75 

0.75 

0.988 

0.836 

1.06 

96 

Orthonol 

MC-1391 

1.00 

0.75 

1.300 

0.870 

1.779 

95 

Orthonol 

Special 

.5916 

.5916 

1.232 

0.905 

0.95 

95 

3.2.3  Fabrication  and  Testing 

As  the  fabrication  of  the  10KW  transformer  proceeded 
and  the  personnel  became  more  proficient  in  the  assembly 
of  the  unique  structures  involved,  several  new  fabrication 
techniques  were  conceived.  The  first  coil,  (Fig  3-8)  was 
constructed  in  a manner  similar  to  the  50  KVA,  400  HZ 
transformer.  It  had  windings  layered  back  to  back  with 
paper  insulation  between.  Each  pair  of  layers  was  then 
separated  from  the  next  pair  by  1/16"  square  wood  spacers. 

The  final  assembly  was  then  vacuum  impregnated  with  varnish. 

It  was  found  that  although  this  coil  occupied  only  slightly 

more  than  half  of  the  core  window,  it  had  wide  variations  < 

in  cooling  passage  dimensions  and  thus  would  tend  to  allow 

non  uniform  cooling.  The  second  coil  wound  (Fig.  3-9)  was 

made  with  individual  single  layers  separated  by  1/16 ' square 

wood  spacers.  A transformer  with  this  coil  was  installed 

in  a case  for  testing  in  Freon  113. 

During  the  period  when  these  coils  were  being  fabricated, 
the  techniques  for  fabrication  of  pie  wound  coils  were  being 
developed.  Difficulties  with  mandrel  alignment,  tensioning, 
precuring,  and  handling  have  all  been  solved.  A complete 
set  of  secondary  windings  have  been  fabricated  as  shown  in 
(Fig.  3-10).  These  are  made  of  #34  AWG  double  coated  wire 
which  can  be  made  to  self  bond  in  a precuring  operation. 

Unfortunately , double  coated  wire  is  not  available  in 
sizes  above  #23  AWG  and  even  this  is  only  available  in 
1,000  lb.  lots.  Therefore,  a new  spoked  mandrel  is  being 
fabricated  which  will  permit  temporary  bonding  of  pie 
wound  coils  with  any  size  wire. 
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Pie  Wound  Secondary  Coils 


In  the  course  of  experimenting  with  pie  winding,  a 
technique  for  multiple  layering  was  developed  which 
permits  much  more  dense  packinq  of  the  wire  than  previous 
methods.  This  technique  will  be  employed  in  the  fabrication 
of  all  pie  wound  coils  in  the  future. 

The  second  layer  wound  transformer  was  installed  in 
a temporary  container  and  successfully  tested  at  full 
rated  current  with  60  HZ  excitation  on  both  primary  and 
secondary  wiiidings  to  cancel  the  flux  in  the  core.  The 
total  loss  was  770  watts  which  is  slightly  above  the 
predicted  loss  of  682  watts.  This  will  be  drastically 
reduced  with  the  new  pie  windings.  The  leakage  inductance 
of  this  transformer  was  found  to  be  25.3  microhenries  which 
is  also  high.  This  will  be  reduced  to  the  desired  level 
of  less  than  4.0  microhenries  by  interleaving  the  pie 
windings . 

3.2.4  Inverter  Cooling  System  - 10KW  System 

The  following  is  a preliminary  cooling  system  design 
for  the  10  KW  inverter  transformer/rectif ier , and  analysis 
of  that  design. 

Input  Data  and  Assumptions 

Transformer 

Copper  Loss  511  Watts 

Core  Loss  59 

Total  Dissipation  570 

Internal  Height  2.30  Inches 

Internal  Depth  2.49 

Internal  Width  3.75  (4.00  with  clearance) 

Rectifiers 

Type:  (4)  Unitrode  687-8  series 

Approx  total  dissipation  = l.lx  13  x 1x8=  114.4 

Overall  Height  2 x 1.14  - 2.28 

Overall  Depth  2 x .75  =1.50 

Overall  Width  3.015 

Note  that  the  height  and  width  package  dimensions  are 
determined  by  the  transformer. 

Total  Power  = 570+114.4  = 684.4,  say  700  watts 
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Package  internal  depth: 


Transformer 
Rectif iers 
Clearances 
Total 


2.49  Inches 

1.50 
.38 

4.38 


Cooling  System:  Use  heat  sink  finned  on  both  sides,  machined 

as  required  to  form  the  top  surface  of  the  T/R  unit.  Cooling 
air  will  be  provided  by  a blower  arranged  for  impingement 
flow  (cooling  air  flow  from  blower  is  normal  to  the  finned 
heat  sink  base  area.). 


For  initial  design  analysis,  use  Thermalloy  extrusion  #6462 
heat  sink  (.250  spaces,  .125  thick  fins,  .655  in  deep  fins). 


Condensing  Side 


Internal  plan  area  of  T/R  unit  = 4.18  x 4.00  inches 
Condenser  area  = no.  fins  x area/fin  + base  area 
No.  fins  in  4 inch  width  = 10 
Area/fin  = 4.38  x 0.655  x 2 = 5.74  sq.in 

Total  area  = 10  x 5.74  + 4 x 4.38  = 74.90  sq  in  0.5  sq.ft. 
Condensing  heat  transfer  coefficient,  hc 


- o 9 4-  3 


o 2 s' 


(22) 


where  F113 

PL  = liquid  density,  lb/cu.ft.  97.69 

PV  = vapor  density,  lb/cu.ft  0.461 

g = gravity  constant,  ft/hr2  4.173  E8 

h = heat  of  vaporization,  BTU/lb.  63.12 

k = thermal  conductivity,  BTUh/(F-ft)  0.038 

mjj  = viscosity,  lb/(ft-hr)  1.646 


then 

K 


2 fe  4 . 8 S 


^ o . 2 *r 


STui/(t1  ft  • F ) (23) 
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For  L = 0.655  inches  = 0.546  ft. 
hc  = 337.08/  a t0*25 

Since  Q = hA  a t, 

(tp— t^)  = Q/hA 

Where 

tp  = fluid  temperature,  F 
tc  = condenser  or  case  temperature,  F 
Q = thermal  load,  BTUh 
A = condenser  area,  sq.  ft. 

tF-tc=  (70 o)  ( 1!/  337.00 

( tF  - tj  = 0.853 

tp  - tt  - 0.84-  F 


(24) 

(25) 


(26) 

(27  ( 
(28) 


External  Forced  Air  Cooling 

Assume  Rotron  "Spartan'1  fan,  type  S,  motor  series  682YS 
Operating  point  = 50  CFM  d>  .175  inches  water 
Forced  air  heat  transfer  coefficient: 

No  fins  on  air  side  (4.687  wide)  = 12 

Effective  area  = 12x. 655x2x4. 687+4. 6872x2/3=  88.3  sq.in. 

Total  flow  area  (two  directions)=  0.25x0.655x12x2  = 3.93  sq  in. 
At  50  CFM,  air  velocity  = V = ( 50) xl44/3 . 93  = 1832  ft/min 
At  assumed  average  air  temp  of  110F, 
air  density  = 0.0696  lb/cu  ft 
air  viscosity  = 0.0465  lb/(ft-hr) 

Equivalent  diameter  = 4x0 . 25x0. 655/ ( 2x (0. 25+0. 655 )) =362  in  = 

.03  ft.  (29) 

Re  = Reynolds  no.  = . 0696xl832x60x . 03/. 0465=4936  (30) 

h = 0.023  k/D  Re°-8Pr°*4  (31) 

= 0.023x0. 0160/0. 03x4936°-8x0.7°*4=  9.61  BTUh/(sq  ft-F)  (32) 


I 
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(33) 


Since  U = hA  (tc  ~ t 


t c ^a.  0 / (7  A) 


(34) 


where 

t = case  temperature,  F 
ta  = air  temperature,  F 
h = heat  transfer  'coefficient , BTUh/(sq  ft-F) 
A = heat  transfer  area,  sq  ft. 


tc  - ta 


(700)  (77)  (3.4'3)/  f(3.6/)(SS.3/  M4\  = 47.7  F (35) 


The  fluid  temperature  will  average  47.7  + 0.9  = 48. 6F 
above  the  ambient  air  temperature,  after  repeatable  cyclic 
conditions  have  been  attained.  Assuming  a maximum  ambient 
temperature  of  100F,  the  fluid  temperature  of  approximately 
150F  will  generate  an  internal  pressure  of  approximately 
15  PS1G. 


A more  detailed  transient  thermal  analysis  will  be 
performed  later  in  the  program  to  determine  the  rate  of 
cyclic  temperature  increase  and  the  peak  temperature 
excursion  over  the  average. 

The  attached  Figure  3-11  indicates  the  general  outline 
of  the  T/R  system.  The  figure  is  a preliminary  concept, 
not  necessarily  to  scale,  and  not  necessarily  reflecting 
the  final  design.  For  example,  the  condensing  fins  appear 
much  longer  than  necessary,  considering  the  low  temperature 
drop  across  the  condensing  surfaces.  It  is  likely  that  these 
fins  will  be  reduced  in  length  in  the  final  design,  reducing 
the  overall  height  somewhat. 

3.3  Dynamic  Analysis  Programs 

3.3.1  Characteristics  of  Magnetic  Cores 

The  directly  measureable  characteristics  of  transformers, 
while  exhibiting  certain  nonlinearities , must  be  altered 
either  by  signal  processing  or  computation  to  yield  information 
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on  the  properties  of  the  feromagnetic  materials.  Each 
manufacturer  has  a preferred  set  of  tests  for  core 
evaluation  depending  on  the  particular  application. 

IEEE  standards  164  and  106  describe  in  some  detail 
the  particular  tests  which  are  deemed  acceptable,  but 
the  degree  to  which  they  are  applied  seems  to  vary  from 
one  manufacturer  to  another.  A list  of  tests  which  have 
been  implemented  at  TTL  is  shown  in  Table  3.  Additional 
tests  for  pulse  and  biased  core  transformers  are  described 
in  references  [1] through[4]. 

TABLE  3 

Magnetic  Device  Tests 

1.  Core  Loss  and  Volt-Ampere  Test  (open  circuit  test) 

2.  Winding  Loss  Test  (Short  circuit  test) 

3.  Sine  Voltage  Tests  (B-H  curve  and  V-I  curve) 

4.  Sine  Current  Tests  (B-H  curve  and  V-I  curve) 

The  specific  materials  tested  at  TTL  have  been  limited 
to  tape  wound  cores  having  lamination  thicknesses  from 
h to  4 mils.  The  materials  used  in  these  tests  were 
Orthonol  (50%  Nickel,  50%  Iron),  Magnesil  (3%  Silicon, 

97%  Iron),  and  Supermendur  (49%  Cobalt,  49%  Iron,  2%  Vanadium). 

A few  samples  of  Permalloy  (79%  Nickel,  17%  Iron,  4%  Molybdenum) 

and  Supermalloy  (78%  Nickel,  11%  Iron,  5%  Molybednum)  were 
tested  along  with  two  ferrite  cores,  but  the  most  extensive 
testing  was  performed  on  the  first  three  materials  mentioned. 

The  purpose  of  these  measurements  was  to  obtain  data 
for  development  of  a mathematical  model  for  transformers 
which  included  the  nonlinear  properties  of  their  cores. 

This,  of  course,  was  not  the  first  attempt  that  had  ever 
been  made  to  model  magnetic  materials  properties.  It  was 
shown  by  Manly  (Ref  6)  that  models  based  on  the  hyperbolic 
tangent,  arctangent,  and  lognormal  functions  fall  short 
of  accurately  representing  the  core  properties.  He  pointed 
out  that  even  models  which  appear  to  closely  fit  measured 
hysteresis  curves  fail  to  match  the  I -V  characteristics. 

Other  models  which  have  been  developed  include  exponential 
series  (Pef  7);  rational  fractions  (Ref  8),  and  least 
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square  curve  fitting  routines  (Ref  9) . Since  none  of 
these  were  easily  adapted  to  the  requirements  of  a 
transformer  analysis  and  design  program  using  minimal 
computing  equipment,  an  effort  was  undertaken  to  develop 
a satisfactory  model. 

Such  a model  has  been  developed  by  TTL  and  its  utility 
in  transformer  analysis  has  been  demonstrated.  The  transition 
from  dynamic  analysis  to  design  has  not  been  made.  So  at 
this  writing,  it  has  not  been  established  that  a transformer 
can  be  synthesized  from  the  TTL  model. 

The  TTL  Model 

The  TTL  model  is  generated  from  three  experimentally 
determined  parameters.  These  are  the  coercive  force,  the 
saturation  flux  density,  and  residual  flux  density.  The 
mathematical  form  on  which  it  is  based  is  y = x/(a+  x ) (36) 

where  x is  a function  of  the  magnetic  field  intensi  ty. 

This  function  has  the  value  zero  at  x = 0 and  approaches 
1 as  x becomes  large  in  either  the  positive  or  negative 
directions.  Its  resemblence  to  a magnetization  curve  was 
one  reason  it  was  chosen  as  the  basis  for  the  model.  The 
mathematical  expressions  for  various  core  properties  are 
shown  below. 

BH  Curves 

The  major  loop  is  modeled  as  shown  in  Equation  (37) 

B - Be,  (H  + K He  )j [He  ~ + | H + K He  (37) 

This  Equation  exactly  matches  experimental  data  at  the 
points  (HcjO),  (0,Br),  (0O,Bs).  The  symbol  k,  has  a value 

of  -1  for  the  lower  or  right  hand  curve  and  +1  for  the 
upper  or  left  hand  curve.  To  represent  a minor  loop,  the 
upper  and  lower  curves  are  shifted  toward  each  other  by 
equal  amounts  such  that  they  intersect  at  the  points 
(j_Hm,  + Bm)  . This  shift  is  represented  by  a term 
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designated  the  displacement  flux  density,  BQ. 


Equations (38) and ( 39) can  be  used  to  generate  a complete 
set  of  B-H  curves  as  shown  in  Fig.  3-12. 


B = + -/)  + |//+AV£|]  -KBk 

8$  *\  + He  ) 


*o  = i 


Hc  (J;  ->)  + \H^HC 

Magnetization  Curve 


8s  (h^-Hc  ) 
Bs 

Bt- 


Hc(f  '0  + 1 


(38) 


(39) 


The  magnetization  curve  Equation  (40)  is  modeled  as 
the  locus  of  the  points  (Hm,Bm)  from  Equations  (38)and  (39) 


Bs  (/-/m  + He) 

Bs  ( Hrr,  - He) 

Hc( % '>)  + Hc  \ 

->)  + \ Hr*-Hc  I 

(40) 


H-VS-B 


The  inverted  form  of  Equation  (38)  is  shown  in 
Equation  (41) 


This  form  is  particularly  useful  for  determination  of  the 
current  from  a known  flux  waveform. 

Core  Loss 

Since  the  model  represents  a continuous  closed  hysteresis 
loop,  the  integral  over  one  complete  cycle  will  yield  an 
expression  for  the  core  energy  loss  per  unit  volume. 


+ Hc  I - \H„  -Hc  I + ch  I — + L - 


C+  I //m  + He 


-4BDH„ 


(42) 
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If  the  parameters  used  for  Bs,Br,  and  Hc  are  for  some 
frequency,  F,  then  the  total  core  loss  will  be  represented 
by  the  product  of  Equation  (42)  and  the  frequency.  For 
dc  data,  Equation  (42)  represents  only  the  hysteresis  loss 
and  the  loss  at  frequencies  other  than  zero  must  include 
a term  for  eddy  current  loss  such  as  shown  in  Equation  (43) . 


Pe  = SffrfT  8^  f <43) 

Permeability 

The  incremental  permeability  may  be  modeled  by 
differentiation  of  Equation  (38).  The  result  is  shown  in 
Equation  (44) . 

§ -')+ <44) 

Since  the  incremental  permeability  approaches  zero  as  the 
magnetic  field  intensity  approaches  infinity,  the  permeability 
of  free  space(/M0  , is  added  to  Equation  (44)  to  model  the 
total  permeability  of  the  material. 

A term  often  used  to  describe  rectanqular  core  materials 
is  the  squareness  ratio,  Br/Bs . This  quantity  approaches  unity 
for  a perfectly  rectangular  core.  The  ratio  of  the 
permeability  at  H = 0 to  that  at  H = He  is 

■S'/Sc  =(/-#f)2  (45> 

From  this  Equation  it  can  be  seen  that  for  a squareness 
ratio  of  0.9  the  slope  of  the  sides  of  the  modeled  loop 
is  one  hundred  times  that  of  the  ends.  Thus  the  model  is 
capable  of  representing  extremely  square  materials,  a 
property  which  is  lacking  in  most  other  core  models. 
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I -V  Characteristics 


Two  limiting  cases  occur  for  a driven  core.  They 
are  voltage  source  and  current  source  conditions,  both 
of  which  are  ideal  cases.  In  Ref.  6,  a sinusoidal  current 
source  given  by 


X - X ^ S/N  uj  t 


(46) 


which  would  force  K to  have  the  form 


K - COS  uj  t 


COS  uj  t 


(47) 


It  may  be  shown  that  the  induced  voltage  for  a core  of 
cross-section,  Ac , and  magnetic  path  length,  , is  given 

by 


y ---  -izY  1481 

The  derivation  of  VI  characteristics  for  a voltage  source 
requires  the  numerical  integration  of  a first  order 
differential  Equation,  and  is  discussed  in  the  section  on 
dynamic  modeling. 

Comparison  of  Model  with  Experimental  Data 

Experimental  and  computed  hysteresis  curves  are  shown 
in  Fig.  3-13  for  Supermendur  at  77°K.  It  appears  that 
the  slope  of  the  sides  of  the  computed  curve  are  less  than 
that  of  the  data.  This  problem  is  due  mainly  to  the  fact 
that  the  value  of  Bs  is  hard  to  determine  since  it  is  not 
possible  to  drive  H to  infinity.  An  alternative  would  be 
to  measure  the  permeability  at  H = Hc,  and  to  compute  the 
saturation  flux  density  from 

Bs  ~ hc  Sf.  ! (/^nc  h/c  -&r  ) 


(49) 
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This  yields  a value  of  Bs  = 1.59T,  whereas  the  value  used 
in  computing  Fig.  3-13  was  1.7T.  The  recomputed  hysteresis 
curve  is  shown  in  Fig.  3-14. 

Data  for  a Magnesil  core  with  the  following  Table  4 
specif ications , are  shown  in  Fig.  3-15. 

TABLE  4 

Magnesil  Core  Type  50086 -2K 
N = 50 

Ac  = 4.3  x 10-6 
& = 8.47  x 10-2 
Br  = 1. 5T 
Hc  = 43.9  at/m 
Bs  = 1.79T 
f = 60  HZ 

The  top  curve  is  hysteresis  and  the  lower  curve  is  the  I -V 
characteristic,  both  were  taken  at  five  different  values 
of  magnetic  field  intensity.  The  computed  BH  curves  and 
I-V  characteristics  are  shown  in  Fig.  3-16.  The  most 
obvious  discrepancy  is  in  the  curvature  of  the  peaks  of 
the  I-V  characteristics.  This  is  thought  to  be  due  to 
the  fact  that  the  source  was  loading  slightly  at  the  peak 
so  that  the  current  was  slightly  non-sinusoidal . This 
theory  is  supported  by  the  analysis  for  the  constant 
voltage  case  which  shown  considerable  rounding  of  these 
curves . 

Examination  of  the  computed  I-V  character istic  shows 
a discrepancy  with  experimental  data  in  that  the  peak 
values  are  shifted  to  the  left  much  more  in  the  data  than 
in  the  computed  curves.  Examination  of  the  B-H  curves 
shown  the  same  phenomenon.  Namely  the  H axis  intercepts, 

Hc,  are  all  much  closer  to  the  Hc  in  the  computed  results 
than  in  the  data.  Experimental  and  computed  values  of  H - 
vs  - Hm  are  shown  in  Fig.  3-19.  From  this  curve  it  is  seen 
that  the  TTL  model  in  its  present  form  deviates  considerably 
from  the  experimental  data.  Residual  magnetisation,  Br,  is 
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Computed  Hysteresis  Curves  for  Supermendur  at  77°  K 
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I- V Characteristics 


Electrical  and  Magnetic  Properties  of  Magnesil 

Figure  3-15 


Ampere  turns  per  meter 
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Hm  - Ampere  turns  per  meter 
Experimental  and  Computed  - vs  - H,^  for  Magnesil 


Fig.  3-17 
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plotted  as  a function  of  maximum  flux  density,  Bm/  in 
Fiq.  3-18.  Here  it  is  seen  that  agreement  between 
theory  and  experiment  is  excellent.  A number  of 
modifications  have  been  attempted  to  correct  the  model 
by  altering  the  Hc  term  in  the  Equation  as  a function 
of  Hm.  The  results  of  two  of  the  more  promising  approaches 
are  shown  in  Fig.  3-19.  It  is  seen  that  the  correction 

has  little  or  no  effect  on  the  curves  while  the  correction 


Ho  = [hc  * j \H * + Hc  ] (51) 

tends  to  reduce  the  H axis  intercept  too  much.  There  is 
some  question  as  to  the  importance  of  this  problem  since 
the  electrical  behavior  is  determined  primarily  by  the 
maximum  flux  density  and  the  instantaneous  permeability. 
These  two  characteristics  are  accurately  represented  by 
the  present  model.  While  high  accuracy  is  desirable  in 
modeling  all  aspects  of  the  material  characteristics  it 
may  not  be  required  in  the  present  program. 

3.3.2  The  Dynamic  Model 


The  basic  nonlinear  transformer  model  is  shown  in 
Fig.  3-20. 


-AyW 


L(-i') 


Fig.  3-20 


Basic  Model  - Core  Only 


The  describing  differential  Equation  is 
lH)  dt  = eU)  -Rl 

where  is  (usually)  a sine  wave  input. 


(52) 

Given  that 
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Modified  BH  Curves  Compared  to  Original  Model 


where 


u>  - 


ZTf 

T ~ Z'W'fo  RAO/ SEC. 


(54) 


we  will  find  it  convenient  to  effect  a change  of  time 
scale  so  that  one  period  of  the  sine  wave  requires  one 
second  in  terms  of  a new  time  variable  T . Let 


ZTTt  = ZTTT 


so  that 


<2?r/0£  = zrrr 


therefore 


r =/0 1 


and 


dt  'Jo 


Equation  (52) becomes 


dL  _ dT 

dt  dt  ~ ^ (r)  ~ Rr 


or 


dZ 

dt 


f(r)  ~ R a. 

/o  £(Z) 
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(55) 


(56) 


(57) 


(58) 


(59) 


(60) 


I J 


n 


Equation  (60)  is  the  basic  Equation  which  was  programed 
for  the  HP  9830  (the  HP  9830  uses  a modified  basic  language) . 


the 


The  nonlinear  model  of  the  core,  Lp£),  is 
equation 


A/  ^ A c q. 

£ 


A.-AcSo.  ^ 

a£  By  1 


defined  by 


(61) 


where,  for  illustrative  purposes,  values  in  Table  5 have 
been  used  as  the  parameters  of  the  model: 


A/  - SO 

Ac  - 4.3*/o  (‘/>nz 
— 8.47  x /O  7 /m 
So  - 4 7T  x /O'7  H//m 

Bs  - /.ST 
Br  = /.ST 
/4c  — 40  Q."t //m 


/4AO//ET/C  CoS£ 
TYPE  50086  -2K 


TABLE  5 

Model  Parameters 
Equation  (60)  can  be  written  as 

^ - £(r>/#  (62) 
dr  -fo'RAC*.) 

and  then  integrated  (graphically)  as  indicated  in  Fig.  3-21. 


Fig.  3-21 

Graphical  Integration  Procedure 
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Equation  (60)  gives  the  slope  of  a solution  curve  in  the 
yi(r)\/sT  plane.  Given  a value  of  current  si  at  T-  % then 
a line  dropped  from  A to  the  curve  ' > has  a length  equal 

to  x.‘  - y . This  defines  the  point  B.  Next,  on  a line 
parallel  to  the  T axis  lay  of  a line  segment  which  is 
equal  in  magnitude  to  ^-L(l)  . This  defines  the  point  C. 
By  construction,  the  line  A C then  has  a slope  equal  to 


'*n 


+ - e(r)/R 

~ I 


(63) 


and  hence,  a "short’  tangential  segment  of  the  line  A C 
constitutes  a solution  to  Equation  (62).  One  then  moves 
to  the  endpoint  of  the  short  tangential  approximation  and 
repeats  the  procedure  to  obtain  the  next  approximation 
to  the  solution  curve.  Given  that  these  tangential 
approximations  are  "short  enough",  one  obtains  an 
approximation  to  the  true  solution  curve  that  can  be  made 
arbitrarily  close. 

The  graphical  procedure  can  be  implemented  as  a 
difference  equation  on  a computer  by  writing 


_ /rr\  ^ AT  ^ AT"  -t 


(64) 


where 


Vo 


sm  , *n  AT)  ~ 


-g  s/>j(nr&Tsrx ) - 

A >2Ac/*o 


h 

R 


N Ac  A?  , £ 


e< 


£ 


-f- 


*(%->) 


dV 


K - ~ ! , \a/H£*  de  ^ G 


(65) 


(66a) 


and 


d jO 


K - / ) dt 


> 0 


(65b) 
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When  an  engineer  graphically  solves  the  problem  using 
Equation  (62)  he  learns  to  adapt  the  "size1  of  the 
tangential  approximation  to  the  steepness  of  the  solution 
curve.  That  is,  for  shallow  slopes  he  is  willing  to  let 
the  tangential  approximation  (the  line  element  CD  in 
Fig.  3-21  be  larger  than  when  the  slope  is  "steep". 

We  can  approximate  the  physical  feel  the  engineer  develops 
(who  is  graphically  solving  the  problem)  by  using  the 
logic  available  on  a digital  computer  to  permit  a variable 
solution  size  that  is  slope  dependent.  For  example,  the 
computer  calculates  a new  point  D based  on  a fixed 
increment.  However,  before  the  point  D is  used  as  a new 
jumping  off  point  in  the  computation  of  the  solution, 
the  slope  of  the  line  CD  is  checked.  If  the  slope 
exceeds  a specified  value,  then  the  AT  interval  is  halved 
and  the  computation  is  repeated  - starting  of  course  from 
the  same  initial  point  A.  This  procedure  is  repeated 
until  the  slope  falls  within  the  constraint  placed  on  it. 
Obviously,  a compromise  must  be  effected-if  we  make  the 
slope  constraint  too  severe  then  the  digital  process  is 
slowed  up  and  one  must  pay  an  excessive  price,  in  terms 
of  computer  time,  for  a degree  of  precision  that  may  not 
be  warranted. 

Given  that  the  current  waveform  must  be  recorded  at 
equal  time  intervals  in  order  to  do  a time  series  analysis 
(for  example,  a spectral  analysis)  then  the  computer  must 
also  be  programmed  to  record  data  only  at  A T second  intervals. 
For  this  reason,  the  procedure  of  halving  the  intervals 
is  a good  one,  since  the  computer  can  remember  the  total 
number  of  times  that  the  interval  was  halved  and  hence 
record  only  the  AT  second  interval  points  for  later  use 
with,  for  example,  a Fast  Fourier  Transform  routine. 

Core  Model  with  Resistive  Load  Ladder  Analysis 

The  model  of  Fig.  3-20  can  be  upgraded  to  that  of 
Fig.  3-22  in  a relatively  easy  manner.  In  Fig.  3-22  let 
the  series  resistance  R become  Rj  and  call  the  pure 
(equivalent)  resistive  load  R2 • 
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Fig.  3-22 


Core  with  Resistive  Load 
Using  Ladder  analysis,  we  can  write 


L(lz')S  + #lJr  L(tz)s  / f0 
1 / 

R, 

A 

U- 

Rz 

SL(Xz)  /£f 

o 

‘ A 

/ 

l(z2)S 

L(z. z)S  / £z 

TABLE  6 

Ladder  Analysis  of  Transformer 
with  Resistive  Load 


From  Table  6 , it  is  seen  that 

-Tz  _ L 


(J^+/)S  L(Tt)+R, 


(6 


or,  in  the  time  domain 


(IH 

d t 

(6 

or 

d A-  2. 

P(t)  Rz  - *'Rz  - Af 
_ ei  ) *.+Ri  R,+Rx  AZ 

(6 

dt 

J-M 

Accounting 

for  the  change  of  time  scale, 

Equation 

(68)  becomes: 


- — . 


(77) 


To  summarize,  the  magnetizing  current  is  found  by 
solving 

<S.(7-J  J$A.  — ^ ^ S- 

■*  a. - _ A?,  r A^ 

for  . The  Load  Current  is: 


~ -^LoiT  (79 

r ^2. 

and  the  total  current  drawn  from  the  source  is . 

Note  that  a "mix"  of  both  frequency  domain  and  time  do- 
main notation  has  been  employed  to  determine  the  governing 
differential  equation.  This  "mix"  of  approaches  is  an  efficient 
analysis  method  and  produces  no  error  since  the  final  result  is 
written  in  time  domain.  In  this  regard,  4.  C-ra)s  has  been 
employed  as  an  operator  and  correct  results  are  obtained  given 
that  jL  CTjJ  s operates  only  on  ta  and  that  the  final  re- 

sults are  written  as  differential  equations  in  the  time  domain. 


Core  Model  With  Resistive  Load  - 


ivalent  Circuit  Analysis 


Fig 3-22  canalso  be  analyzed  by  the  use  of  equivalent  cir- 
cuits. First,  convert  the  voltage  source  to  a Norton  equiva- 
lent current  source  (Fig. 3-23). 
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Fig. 3 -2 3Analyses  Using  Equivalent  Circuits 


Next  find  the  parallel  equivalent  of  and  /TA_(Fig.3-23h)and  th 
convert  back  to  an  equivalent  voltage  source  (Fig.  4c) . From 
ig.3  23C  one  obtains  directly  the  differential  equation  des- 
cribing the  magnetizing  current  as 

e*-  --  ..  /ZS.'?*., (81) 

which  verifys  Equation  (68 ) . 

Additional  Dynamics  - Series  Leakage  Reactance 


The  next  step  is  to  modify  the  circuit  to  account  for 
an  equivalent  series  leakage  reactance.  (Refer  to  Fig. 3-24  where 
has  been  added.) 


Fig  .3  -24  Leakage  Reactance  Added 
Ladder  analysis  can  now  be  applied,  given  that  the  proper 
interpretation  of  , as  an  operator  on  alone, 

is  observed. 
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Table  7 Ladder  Analysis  - Series  Leakage  Reactance 

From  Table  7 .the  ratio  of  (the  magnetizing  current)  to  £~u 

is 


Jf  = 


^sj  '^‘■a. 


Let 


‘Zy  ~ 


Sa>^- 


5r '" 


^F*j  ^ . 


Iss?/" 


1 


[«^] 


The  equivalent  circuit  is  given  in  Fig.  3_25. 
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Fig. 3-25  Equivalent  Circuit 
Let  -T”'  [^1  = 

^-/  [^]  * ev^ 


(8  6) 
(87) 


so  that  the  differential  equation,  describing  the  magnetizing 
current, becomes 

jLs^)  £yCc-4> -*4.00*1*  (ae) 


fr-xJ-fxC*)***  (89.) 

The  precise  form  of  Equation  (8  9)  depends  on  the  assumed  values 
of  the  parameters.  For  example,  as  ->■  ^cij  an<*  ' 

On  the  other  hand,  when 


or 


- gy.<»V$y 


*/ 


(90) 


For  this. important . case.  Equation  ( reduces  to 

d*  ■£('*■*■)  (9  2) 

and  where  «^>^is  dependent  on  the  form  of  the  assumed  voltage. 
Given  that  then 

L **/  J 

= ^ , / (94) 


/9-5 

= 3*-rU>A  S*+Cu 


Therefore, 


— >9  '«*•  -t-  4*>  £ »*  <^e.  "X* 

/9  = -C  <=  ^4/  ^ ^ *) 

^ =.  vs  f j f 


In  general,  will  still  have  the  form  given  in 

Equation  (96)  when  both  and  ^ are  retained.  However, 

the  coefficients  of  the  partial  fraction  expansion  will  be 
different.  Specifically, 


— C = 


V.  ^ 


: _£A£±f, 


- 

&^jS=j!  J 


Note  that  the  increase  in  the  "order"  of  the  equivalent  voltage 
source,  due  to  the  introduction  of  the  energy  storage  element 
does  not  particularly  complicate  the  problem.  That  is, 
the  effect  of  L,  on  eptfcan  be  treated  in  an  analytical 
fashion  - there  is  no  need  to  "simulate"  the  complicating 
effects  of  L,  as  far  as  the  equivalent  source  is  concerned. 

The  effect  of  on  the  equivalent  impedance  is, 

however  more  involved  and  does  theoretically  require  that  we 
increase  the  order  of  the  simulation  model.  Referring  to 
Equation  (84 ) , observe  that  the  equivalent  impedance  approaches 
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a constant  value,  for  all  frequencies,  only  for  two  cases: 
a)  -f-y » b)  . For  the  intermediate  range  where 

^ it  becomes  necessary  to  find  the  exact  "impulse  re- 

sponse" of  2(s>  • To  find  out  what  this  impulse  response  is, 
write  Equation  (84)  in  rational  form  (i.e.,  the  ratio  of 
polynomials  with  the  numerator  polynomial  of  one  order  less 
than  the  denominator) : 


, _ -e*-  [JzisA . ^ 

J 4., 


pod 


= ^ 


S+ 


s&~  % c 


(102) 


(103) 


Using  Equation  (103^and  Equation  (96)  , Equation  l89) 

becon,esv,^  -fiksy-lzA 

**  (104) 


Equation  ( 104)simplifies  only  slightly: 

'V'i7 


(105) 


(106) 


so  that 


.2  - u/ 

<<?<«■  z ^ z. 


(107) 


Thus  Equation  (105) can  be  replaced  by  the  system  of  equations 


<Z  ■<  j . ^ &xi ■ ^ ^ 

' ^7~ 

v.  * -s  - #,)  u/ 

: 4 4 k~ 


r'ftC* 


nos) 


(109) 


This  system  of  equations  has  been  programmed  for  the  HP9830 
and  works  well  when  and  Z has  a reasonable  numerical 
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value  - on  the  order  of 


order  of  /o'8  henrys  the 
Equation  (96)  is  recommended  only  when  has  a significant  value. 
This  condition  can  be  checked  quite  nicely  using  Equation  (84). 
Note  that  Equation  (84 ) gives  the  form 


/0  henrys.  However,  for  A on  the 
HP9830  underflows.  Hence  the  use  of 


r 


• 

~Z7 


s + / 


s+. 


or 


7' 


^ : 

r-S-y-/ 

1 

A". 

( 110) 


( ill) 


u>. 


We  may  compare  the  break  frequencies 
frequency  of  the  input  sine  wave  . 

we  may  as  well  set 

^7~ 


and  against  the 
If  then 


( 112) 


! 


In  this  event.  Equation  (89)  becomes 

~zr 274; <113) 


Therefore,  the  program  devised  for  Figure (Equation  (00) 
can  be  used  by  letting 


e&- 


-*Z  * 


(114) 


where 


/? 

A*  -c  = 


i/  r 


(115) 

(116) 
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A More  Detailed  Model 


Consider  Fig  .3-26  which  makes  provision  for  capacitive 
reactances  as  well  as  inductive  losses.  (We  draw  a frequency 
domain  circuit  diagram.) 
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Fig.3-26A  More  Comprehensive  Equivalent  Circuit 

This  circuit  can  be  reduced  to  the  nonlinear  first  order 
form  via  the  successive  transformations  given  in  Fig.  3-27 
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Fig. 3-27  Reduction  of  Fig To  First  Order  Form  of  Fig.  3 20 
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Fig .3-27  (con*  t)  Reduction  of  Fig.3-26to  First  Order  Form  of  Fig.  3 
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gives  the  first  order  form 

_ £ *«<*) ~ / (*'  I'JJ-f-r) 

**  2 

(124) 

Substituting  the  equivalent  impedances  into  Equations 
(12  3)  gives  explicit  expressions  for  .Z^^and  • 

( 122)and 
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One  may  evaluate  quation  (125)  and  (l26),  for  given  para- 
meter values,  and  check  for  resonances,  the  relative  importance 
of  each  term  in  and  £Tej.  , etc.  For  example,  the  importance 

of  the  various  break  frequencies  encountered  in  and 

can  be  assessed  using  Bode  plots,  on  the  basis  of  their 
roltivance  to  the  frequency  of  the  input  sinusoidal.  For 
example,  if  Equation  125  were  to  yield  a break  frequency  at 
5000  /7Z  we  would  choose  to  ignore  it  if  the  input  sinusoid 

were  at  60  //Z.  . However,  a break  frequency  at  5000  would 
have  to  be  retained  in  the  model  if  the  input  frequency  were 
1000  //Z.  These  points  are  demonstrated  with  an  example.  Let 
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Using  these  numbers,  it  can  be  shown  that  Z ^ 


and  Ztp  can  be 


factored  explicitly  as 
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Since  the  input  frequency  is  377  rad/sec,  the  break  fre- 
quencies of  Equations  (130)  and  (131)  are  of  no  concern  and  we  may 
as  well  set 


e*. 

7 *, 

and 
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Load  current  and  magnetizing  current  are  plotted  in  Figs. 3-28 
and  3-29  for  two  different  values  of  input  voltage  (V=3  and 
V=5).  In  addition,  the  spectral  content  of  the  load  current, 
for  the  V=5  case,  is  examined  in  Fig.  3-30 

In  both  cases  (V=3,  V=5) , the  following  parameter  values 
were  used  for  the  nonlinear  model: 
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TABLE  8 


Model  Parameters,  Illustrative  Example 
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Note  the  large  "inrush"  component  in  the  magnetizing 
current  (for  the  3 volt  case)  that  dies  out  to  an  extremely 
small  value  after  the  first  cycle.  This  is  not  the  Cuse  for 
the  Vc  =5  volt  input  - there  it  is  seen  that  the  "inrush" 
current  component,  although  significant,  does  not  dominate  in 
the  plot  of  magnetizing  current  vs.  t'  . 


Note  that  the  frequency  dependent  terms  in  Equations  (125) 
and  ( 126  exhibit  a pole-zero  cancellation  trend  which  seems  to 
indicate  that  only  a few  of  the  time  constants  should  be  re- 
tained in  a much  more  simplified  model  of  the  transformer. 

That  is,  the  model  of  in  Equation  (125) was  a fourth 

order  polynomial  which  reduced  to  a second  order  polynomial 
over  a first  order  polynomial,  for  this  particular  set  of  para- 
meter values,  we  decided  to  ignore  the  additional  dynamics 
introduced  by  this  lead-lag-lead  factor  only  because  the  break 
frequencies  were  sufficiently  high  with  respect  co  the  60  'Aa 
input  sine  wave.  (377  rad/sec).  However , these  additional  dy- 
namics should  be  taken  into  account  for  input  waveforms  that 
have  frequencies  greater  than  about  600  rad/sec. 
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To  conclude,  it  is  apparent  that  the  importance  of  the 
frequency  dependent  terms  in  Equations  (125) and  (126) should  be 
evaluated,  for  each  individual  case  in  order  to  gain  insight 
into  the  character  of  Zn<j  and  . For  example,  using  the 

numbers  given  in  the  illustrative  example,  it  was  seen  that  all 
the  roots  of  the  polynomials  were  first  order  terms.  However, 


3.4 


200KW  Transformer 


The  original  plan  was  for  the  200KW  system  to 
be  a scaled-up  version  of  the  10KW  system.  At  the 
time,  it  was  not  known  what  the  effect  of  nonlinear 
characteristics  would  be  on  this  intended  scaling. 

In  fact  one  of  the  program  objectives  was  to  determine 
the  extent  and  predictability  of  these  nonlinearities. 

As  the  work  progressed,  it  became  evident  that  direct 
scaling  was  not  possible.  This  being  the  case,  some 
preliminary  200KW  designs  were  run  with  resulting 
system  specific  weights  of  less  than  0.1  lb/KVA.  The 
10K  W test  results  will  contribute  to  the  establishment 
of  coil  separation  and  interleaving  requirements  of  the 
200  KW  unit.  Detailed  design  of  this  unit  will  follow 
and  final  10KW  transformer/rectifier  fabrication  and 
testing . 

4.0  Conclusions 

The  program  is  proceeding  on  schedule.  Valuable 
information  has  been  obtained  from  the  10KW  fabrication 
and  testing  which  has  permitted  the  elimination  of 
several  inaccuracies  in  the  computer  programs,  and  the 
addition  of  certain  new  features  such  as  pie'’  winding 
design  capability.  New  techniques  for  the  fabrication 
of  lightweight  'pie "windings  have  been  developed  which 
show  great  promise.  Several  prototype  10KW  units  have 
been  completed  and  subjected  to  preliminary  testing.  A 
10KW  water  cooled  T/R  unit  will  be  completed  in  August  1976 
as  scheduled. 

Several  improvements  were  made  in  the  magnetic  core 
model  which  permit  the  accurate  simulation  of  hysteresis 
and  I -V  curves  for  most  core  materials.  A procedure  was 
developed  for  determining  the  complete  dynamic  behavior 
of  any  transformer  which  can  be  represented  by  lumped 
parameters  and  the  nonlinear  core  model.  This  procedure 
yields  both  time  dependent  characteristics  and  harmonic 
content.  It  requires  that  all  of  the  transformer 
characteristics  be  specified. 

The  200  KW  transformer  design  should  proceed  without 
difficulty  after  the  10KW  unit  is  complete. 
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